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ABSTRACT 


Several  heuristics  and  algorithms  have  been  developed  to  find  minimal  sum-of- 
products  expressions  in  binary  logic.  Most  of  them  use  prime  implicants  during  min¬ 
imization  process. 

An  efficient  search  strategy  has  been  developed  for  finding  minimal  sum-of-products 
expressions  for  multiple-valued  logic  (MVL)  functions  by  using  the  constrained 
implicants  set  concept.  The  search  space  can  be  considerably  reduced  over  the  only 
other  known  exact  minimization  technique  and  exhaustive  search. 

The  primary  goals  of  this  research  are  to:  (1)  examine  whether  the  constrained 
implicant  set  concept  can  be  efficiently  used  in  binary  logic,  and;  (2)  develop  a  heuristic 
called  the  constrained  implicant  set  heuristic  (CISH).  The  general  idea  of  the  CISH  is 
to  select  the  minterm  with  the  least  implicant  cover  size  and  find  the  implicant  with  the 
largest  minterm  coverage  that  covers  a  selected  minterm. 

In  this  research,  the  implementation  of  the  CISH  is  presented,  the  performance 
analysis  of  the  CISH  is  shown  by  comparing  with  other  heuristics  (Maximum  Implicant 
Heuristic,  Espresso  II)  with  respect  to  the  average  number  of  the  product  terms,  the 
average  computation  time,  and  the  average  memory  usage. 
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I.  INTRODUCTION 


A.  MOTIVATION 

The  primary  goal  of  good  design  in  binary'  logic  is  to  find  a  realization  of  a  given 
function  at  minimal  cost.  There  have  been  several  different  costs  in  binary  logic  design. 
The  number  of  the  gates  is  a  commonly  used  cost.  Another  cost  is  the  number  of 
product  terms  used  in  the  sum-of-products  expression  of  the  given  function.  This  crite¬ 
rion  has  become  especially  important  in  recent  years  with  the  introduction  of  PAL's 
(programmable  array  logic)  and  I’LA's  (programmable  logic  arrays).  There  has  been  a 
corresponding  increase  in  interest  in  algorithms  for  finding  the  fewest  number  of  product 
terms  needed  to  realize  the  given  function. 

An  efficient  search  strategy  has  been  developed  for  finding  a  minimal  sum-of- 
products  expression  by  using  the  constrained  implicant  set  concept  in  multiple-valued 
logic.  The  search  space  can  be  considerably  reduced  over  the  only  known  exact  min¬ 
imization  techniques  by  using  the  constrained  implicant  sc-t  concept  [Ref.  1  j.  A  primary 
motivation  of  this  thesis  is  to  see  if  the  constrained  implicant  set  concept  can.  be  ex¬ 
tended  to  binary  logic, 

A  new  heuristic  for  binary  functions  is  described  which  is  called  the  constrained 
miplicam  set  heuristic  (C1S11).  An  analysis  of  the  performance  of  the  CISH  has  been 
done  by  comparing  the  results  of  the  average  number  of  sum-ol-products  term,  compu¬ 
tation  time,  and  memory  usage  with  two  existing  algorithms  1)  Maximum  Implicant 
Heuristic  (Mill)  and  ?)  Espresso  11. 


B.  NATURE  OF  THE  PROBLEM 


The  goal  of  logic  minimization  is  to  find  a  minimal  sum-of-products  expression  of 
a  binary  function.  This  problem  has  received  considerable  attention  for  some  time. 
Early  methods  such  as  Quine-McCluskey  (Ref.  2},  and  iterated  consensus  [Ref.  3]  begins 
the  minimization  by  finding  all  prime  implicants  of  the  function.  The  finding  of  prime 
implicants  is  used  by  most  heuristics  minimization  methods. 

The  interest  in  heuristic  methods  is  due  to  the  large  computation  times  required  by 
exact  minimizations  algorithms.  I  or  example,  a  10  variable  binary'  function  can  have 
as  many  as  5904  prime  implicants  while  a  function  with  20  variables  can  have  as  many 
as  174,339,220  prime  implicants.  'I  he  relationship  between  the  number  of  the  prime 
implicants  and  number  of  the  variables  (/;)  has  been  shown  as  3 "/«  in  the  worst  case, 
(Refs.  4,  5:  p.  49] 

It  is  possible  to  find  an  exact  minimal  sum-of-products  expression  for  a  logic  func¬ 
tion  with  a  small  number  of  variables  or  simple  functions  with  larger  number  of  variable. 
When  the  number  of  the  variable  increases  or  the  function  becomes  more  complex,  then 
mere  computation  time  is  needed  to  extract  the  exact  minimal  solution.  Sometimes 
finding  nearly  minimal  solution  of  a  function  in  shorter  computation  time  has  more  im¬ 
portance  and  advantage  than  finding  an  exact  minimal  representation  in  very  long 
computation  time  [Ref,  o]. 

The  exact  minimal  sum-of-products  expression  can  be  solved  by  enumeration.  In 
■this  method,  all  possible  solutions  aie  tried  and  then  the  one  having  the  fewest  number 
of  sum-of-products  is  chosen.  This  method  needs  very  long  computation  time  due  to 
large  number  of  possibilities  of  potential  solutions.  Besides,  it  has  been  shown  that  ex¬ 
tracting  minimal  sum-of-pioducts  solutions  from  a  complete  set  of  prime  implicants  is 


an  NP-hard  problem.  The  best  known  algorithms  for  such  prob’ems  require  exponential 
time  [Refs.  1,7:  p.  246). 

C.  THESIS  OUTLINE 

Notation  and  definitions  are  given  in  Chapter  II.  The  constrained  implicant  set 
heuristic  is  introduced  in  Chapter  III.  Chapter  IV  and  Chapter  V  discuss  the 
comparision  results  and  performance  with  Maximum  Implicant  Heuristic  and  Espresso 
II. 
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II.  BACKGROUND  AND  DEFINITIONS 
A.  DEFINITIONS  IN  BINARY  I '  GIC 


In  this  section,  we  briefly  summarize  fundamental  definitions  used  in  binary  logic 
(Refs.  2,3]. 

Definition  1: 

A  literal  is  a  variable  or  the  complement  of  a  variable.  Examples:  vv,  x,  vv,  x 

Definition  2: 

A  product  term  is  a  single  literal  or  a  Boolean  product  of  the  literals.  Examples: 
ivy;,  xyl,  vv; 

Definition  3: 

A  minterm  is  a  product  term  where  a  literal  of  each  variable  appears  exactly  once. 
Examples:  vv.vvz,  wxy:,; 

Definition  4: 

A  logic  function  /•’(*„; x2,  covers  a  logic  function  P(  jc„  x>,  ... ,  x„)  if  for  every 

input  combination  such  that  P-  1,  then  F=  1. 

Definition  5: 

Let  1  be  a  product  term  of  function  f.  If  a /is  nonzero  for  all  minterms  covered  by 
I,  then  I  is  an  implicant  of  the  function.  Examples:'  vv;,]7;  •• 
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Definition  6: 

Let  I  be  an  implicant  of  the  function/  I  can  be  said  prime  implicant  of  the  function, 
if  it  is  not  covered  by  any  other  imphcants  of  the  function  / 

Definition  7: 

A  minterm  is  said  to  be  distinguished- 1 -cell  of/  if  this  minterm  is  covered  by  only 
one  implicant  of  the  / 

Definition  8: 

A  prime  implicant  I  is  said  to  be  an  essential  prime  implicant  if  it  covers  a 
distinguished- 1  -cell  of /. 

Definition  9: 

Let  A™,  Aon ,  Apc  be  sets  of  assignments  of  values  to  variables  of  function / respect 

to  minterm  values 

•  A0,v  is  the  set  of  minterms  of / 

•  A0fF  is  the  set  of  assignments  of  values  to  the  variables  such  that /is  0. 

•  ADC  is  the  set  of  assignments  of  values  to  variables  such  that  /is  don't  care. 

Definition  10: 

A  function  /  can  be  considered  to  have  a  cycle  if  it  has  more  than  one  minimal 
sum-of-products  expression.  Example:  see  the  function  in  Figure  6  on  page  17.; 

Definition  11: 

Let  v  and  /?  be  n unterms  such  that  complementing  one  literal  in  a  jields  /?.  Minterm 
a  and  /»’  aie  called  Direct  .Xcighbui s. 
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B.  BACKGROUND  IN  MULTIPLE-VALUED  LOGIC 


The  constrained  implicant  set  concept  was  originally  developed  from  the  work  of 
finding  absolute  minimization  in  multiple- valued  logic  by  Jon  T.  Butler  and  P.  Tirumalai 
[Ref.  1],  Absolute  minimization  tries  to  find  the  absolute  minimal  realization  of  a 
function  by  doing  an  exhaustive  search  of  all  possible  solutions.:  An  algorithm  for  ab¬ 
solute  minimization  is  introduced  in  Appendix  A.: 

As  stated  before,  absolute  minimization  needs  considerable  computation  time.  The 
search  space  is  also  very  large  in  absolute  minimization.  This  space  can  be  made  smaller 
by  applying  a  limitation  rule:  constrain  some  implicants  and  establish  the  constrained 
implicant  set  to  find  minimal  solution.  The  constrained  implicant  set  concept  signif¬ 
icantly  reduces  search  space  as  well  computation  time  to  get  exact  minimal  solution  in 
multiple-valued  logic  [Refs.  1,8]. 

Definition  12  : 

K  (a)  is  a  constrained  implicant  set  of  minterm  a  for  function /,  if 

R(a)  -  U  I 


where  I  is  an  implicant  of /.  [Ref.  1] 

Lemma  1  : 

If  R(a)  is  a  constrained  implicant  set,  then  even-  possible  sum-of-products  ex 
pression  for/has  to  contain  at  least  one  implicant  in  IR(a)  [Ref,  1]. 

Definition  13: 

R{o')  is  a  minimal  constrained  implicant  set  of  a  function  if  and  only  if 
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0  <  |  U?(a)  |  <,  |  IR(/?)  | 


for  all  other  minterms  /?  of f  where  i(a)  is  a  constrained  implicant  set  [Ref.  1], 

The  search  space  can  be  represented  as  a  tree  where  each  node  represents  a  function 
and  each  edge  corresponds  to  implicant  of  the  upper  level  function.  The  root  node  is  the 
given  function  to  be  minimized.  The  functions  at  the  next  level  down  from  the  root  node 
can  be  obtained  by  subtracting  an  implicant  from  the  root  node.  Subtracting  an 
implicant  from  a  function  corresponds  to  setting  1  's  or  don't  cares  in  the  function  cov¬ 
ered  by  the  implicant  to  don't  cares. 

If  the  root  function  has  £  implicants,  the  root  function  has  £  branches  or  subf¬ 
unctions.  Further,  each  subfunction  has  a  maximum  of  £-1  descendents.  It  can  be  seen 
that  when  the  root  function  has  many  implicants,  tire  search  space  is  large.  This  situ¬ 
ation  can  produce  large  computation  time  needed  to  find  the  minimal  solution  of  the 
function.  That  is,  the  solution  is  to  try  each  possible  path  on  the  search  space  where  the 
shoircst  path  (having  the  fewest  number  of  the  implicant)  is  chosen  as  the  minimal  ex¬ 
pression  of  the  function  [Refs,  1,8]. 

By  using  Lemma  I,  the  search  space  can  be  made  smaller.  At  least  one  implicant 
from  the  constrained  implicant  set  has  to  be  in  the  minimal  sum-of-products  expression 
of  the  function.  When  R(a)  is  chosen  as  small  as  possible,  there  arc  fewer  choices  than 
for  larger  size  R(a) ,  This  decreases  the  computation  time  because  there  are  fewer  paths 
to  be  examined,  [Ref.  8] 

Example  1: 

This  example  illustrates  the  search  space  of  a  specific  function  and  the  finding  of 
exact  minimal  solution  of  a  /  in  search  space  by  using  absolute  minimization  algorithm 
and  constrained  implicant  set  concept. 
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A  five  variable  function  with  sixteen  min  terms  is  chosen.  The  function  and  exact 
minimal  solution  is  shown  in  Figure  1. 

Let 


J{vtw,X)y,z)  =  y m{  1,5,  6,  7. 11,  12, 13,  15,  16,  17,  18,22,24,26,27,  28) 

immmJ 

The  numbers  enclosed  in  parenthesis  on  the  right  hand  side  correspond  to  minterms  in 
the  binary  representation.  For  example,  minterm  a  =15  corresponds  to  assignment 
01111.  Specific  piimc  implicants  in  an  exact  minimal  solution  arc  represented  by  a 
capital  letter  as  seen  in  Figuic  1. 

'lhc  search  spaces  for  the  gisen  function  aic  illustiatcd  in  Figuic  2  on  page  10  and 
Figure  3  on  page  11,  Because  of  the  dilliculu  of  showing  all  possible  search  paths,  only 
one  search  path  is  shown  as  example  (that  which  gives  the  exact  minimal  solution). 
Each  node  function  is  indicated  b>  a  solid  dot  in  search  space,  'lhese  node  functions  can 


8 


fv"'  ‘  *-  ’  .  -  J 


be  found  by  subtracting  an  implicant  of  the  function  from  the  root  function  as  described 
on  page  7. 

The  node  functions  on  the  search  path  chosen  as  an  example  are  represented  by 
letters  NF,  These  functions  at  the  next  level  down  from  the  root  node  can  be  obtained 
by  subtracting  the  implicant  (that  shown  by  capital  letter  in  Figure  1)  from  the  root 
node.  These  node  functions  are  stated  at  the  bottom  of  the  Figure  2  and  Figure  3. 

The  search  space  to  find  exact  minimal  solution  /is  shown  in  Figure  2.  The  given 
function / has  fourteen  prime  Implicants,  and  thus  it  has  fourteen  branches  from  the  root 
as  described  in  absolute  minimization.  The  search  space  is  very  large,  wide  and  deep, 
because  there  are  many  node  functions  and  subbranches  in  the  search  space  of/  (U\, 
2,162,160  possible  branches).:  All  possible  search  paths  must  be  investigated  at  each  in¬ 
dividual  node  to  find  the  exact  minimal  solution. 

On  the  other  hand,  the  constrained  implicant  set  concept  can  be  applied  to  /.  A 
search  space  has  been  created  for  the  same  function  as  shown  in  Figure  3  on  page  11. 
This  search  space  has  fewer  branches  (i.e.,  8  )  and  fewer  node  functions  (i.e,,  384).  De¬ 
finitions  and  rules  to  create  a  search  space  and  finding  the  minimal  solution  is  discussed 
in  Section  C  and  Chapter  III, 

As  a  results  of  the  comparision  of  the  two  search  spaces  in  Figure  2  and  Figure  3, 
we  sec  that  although  the  same  implicant  are  chosen,  the  absolute  minimization  algo¬ 
rithm  has  many  more  branches  at  each  node  function.  Applying  the  constrained 
implicant  set  concept  to  the  minimization  of  function  reduces  the  number  of  branches 
and  node  functions  to  be  examined.  Therefore,  the  program  needs  less  time  to  find  the 
minimal  solution. 
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C.  DEFINITIONS  USED  IN  C1SH 


In  this  section,  the  fundamental  measures  used  in  CISII  is  defined. 

Definition  14: 

Let  /CS(a)  denote  implicant  cover  size  of  minterm  a.  The  implicant  cover  size  is  the 
number  of  the  implicants  that  cover  minterm  a. 

When  /CS(a)  is  calculated  for  minterm  «,  the  minterm  a  is  not  counted  as  an 
implicant  in  CISH.  On  the  other  hand,  there  is  an  exception.  That  is,  if  a  minterm  does 
not  have  any  direct  neighbor,  it  is  counted  as  an  implicant  in  calculating  ICS( a). 
Therefore,  the  minterm  with  no  direct  neighbor  has  same  /CS(a)  (i.e.,  1)  as  the  minterm 
with  only  one  direct  neighbor. 

Definition  15: 

Let  A/C(l)  denote  the  minterm  coverage  of  implicant  I,  The  minterm  coverage  is  the 
number  of  the  minterms  that  are  covered  by  I.: 

MC( I)  is  used  to  determine  the  importance  or  the  cost  of  I  during  the  minimization 
process.  The  cost  of  the  implicant  is  introduced  in  the  minimization  as  the  number  of 
minterms  covered  by  1,  The  implicant  with  the  highest  A/C(I)  is  chosen  in  the  minimi¬ 
zation  process. 

In  CISH,  each  minterm  covered  by  a  selected  implicant  I  turns  to  a  don't  care  term. 
Each  new  don't  care  term  covered  by  I  is  subtracted  from  A/C(I).  So  the  don't  care  term 
affects  the  MC( I).  Therefore,  the  MC( I)  changes  during  the  minimization  process  in 
CISH.; 
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D.  MEASURES  USED  IN  CISH 


1.  Clustering  Factor 

7CS(a)  is  a  measure  of  the  degree  to  which  other  minterms  cluster  around  «. 
It  shows  how  many  minterms  with  which  a  minterm  a  can  combine.  The  lower  7CS(a) 
is,  the  fewer  combinations  exist. 

•  j 

Example  2: 

To  illustrate  these  definitions,  consider  the  following  four  variable  function  with 
ten  minterms.  Let 


/(hvw)  =  ^T/»(  0.  3,  4,  7,  8, 9,  10, 1 1,  14.  15 ) 

The  function /is  illustrated  in  Figure  4  on  page  14. 

The  lCS(o c)'s  for  each  minterm  in  Figure  4  are  shown  in  the  corresponding  up¬ 
per  left  corner.  For  example,  the  minterm  «  =  0  has  7C5(a)=  2.  Two  implicants  (I,,  I2) 
cover  this  minterm,  where  1,=  vvyz  and  I2  =*  xyz  .  In  general,  minterms  in  the  center 
of  a  cluster  have  a  higher  7CS(a).  For  example,  minterm  a=  11  is  in  the  middle  of  a 
cluster  of  l's  and  it  has  a  high  7CS(a)  (7),  while  minterm  a  =  4  is  remote  and  has  a  low 
lCS(cf.)  (i.e.,  1). 

2.  Use  of  lmplicants 

As  an  experimental  result,  using  all  implicants  rather  than  prime  implicants  of 
function  /  in  calculating  7C$(a)  provides  more  information  about  minterms  and 
implicants.  We  believe  that  only  using  the  prime  implicants  set  of / may  miss  some  in¬ 
formation  about  clustering  of  minterm  a  with  neighboring  minterms.  Therefore,  all 
implicants  of  given  function  are  counted  in  ICSio). 
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Figure  4.  An  Example  on  Calculation  of  1CS{ a) 


Example  3: 

This  example  illustrates  the  advantage  of  the  using  implicates  versus  ptime 
implicates  to  calculate  ICS(a),  The  same  function  in  Example  1  is  used.  The  minterm 
located  at  00101  (i.c.,  a=5)  is  chosen  as  sample  minterm  to  show  changing  of  1CS( «), 
The  ICS(  a)  of  each  minterm  that  is  counted  by  vising  both  prime  implicates  (Figure  5. a) 
and  all  implicates  (Figure  5.b)  arc  shown  in  page  15- 

In  Figure  5. a,  only  prime  implicates  ate  considered  to  calculate  /CS(a).  All 
minterms  have  same  ICS{ a)  (i.c.,  2).  It  shows  that  all  mintenns  have  the  same  clustering 
with  their  neighbors.;  Thus,  all  ul  them  all  equally  likely  candidates  as  the  starting  point 
for  the  minimization. 

On  the  other  hand,  if  all  possible  implicates  are  counted  in  ICS(v)  we  have  the 
the  situation  shown  in  I'iguie  5.b.  In  that  ease,  the  /CS(p)  of  each  minterm  is  not  the 
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same  as  in  Figure  5.a.  The  ICS(<x )  shows  more  accurately  the  clustering  and  coupling 
strength  among  the  minterms.  Besides,  the  importance  of  the  minterms  to  be  starting 
point  for  minimization  is  significantly  changed  (i.e.,  7CS(a  =  5)=4).  The  minimization 
can  start  any  one  of  the  eight  minterms  where  their  /CS(a)  =  2.  These  minterms  are 
a*  1,  6,  11,  12,  17,  22,  27,  28. 

There  is  a  special  case  that  either  use  of  all  implicants  or  prime  implicants  can¬ 
not  provide  any  advantage  over  other  set  in  calculating  ICS{a).  This  case  occurs  when 
the  set  of  all  implicants  and  is  identical  to  the  set  of  all  prime  implicants;  that  is,  when 
each  implicant  of  the  function  is  a  prime  implicant. 

Example  4: 

This  example  illustrates  the  special  case  in  calculation  of  ICS(a).  A  four  vari¬ 
able  function  with  eight  minterms  is  used  as  an  example.  It  is  shown  in  Figure  6  on 
page  1 7.  ICS(tx)  of  each  minterm  in  Figure  6  is  shown  in  the  corresponding  upper  right 
corner. 

Let 


</(vv,x4y,z)  =  £m(  0,  2, 4,  5, 10, 1 1 , 1 3, 1 5 ) ,: 

It  can  be  seen  that  the  ICS(a)  is  the  same  for  each  minterm  whether  all  implicants  or 
prime  implicants  are  counted  in  lCS{a).  Thus,  there  is  no  advantage  in  using  prime 
implicant  or  all  implicant  in  calculating  7CS(a), 

3.  Effects  of  Don't  Care  Terms 

From  the  definition  of  ICS(u),  all  implicants  should  be  considered  in  counting 
/CS(a).  If  it  is  desired  to  find  the  /CS(a)  during  the  minimization  process,  it  is  necessary 
to  consider  all  don't  care  terms  as  1  terms..  Changing  the  don't  care  terms  to  1  terms  in 
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Figure  6.  Special  Case  in  Calculation  of  /CS(«) 

substeps  of  tltc  minimization  gives  the  original  root  function  (that  we  count  1CS{ a)  at 
the  very  beginning).  Therefore,  /CS(a)  of  each  minterm  will  remain  unchanged  through 
the  minimization  process.  Thus,  it  is  enough  to  calculate  the  1CS{ a)  of  each  minterm 
once  at  the  beginning  of  the  minimization.  It  means  that  the  don't  care  terms  that  in¬ 
troduced  during  the  minimization  process  don't  affect  the  /CS(k)  of  each  minterm. 

i\/C(l)  is  allected  by  don't  care  terms  introduced  during  the  minimization 
process.:  The  unpoituncc  of  the  implicant  is  inversely  proportional  to  the  number  of 
minterms  changed  to  don't  care  at  each  node  function  in  the  minimization,  l-ach 
minterm  changed  to  don't  care  and  covered  by  1  1  educes  the  importance  (increases  the 
cost)  of  the  1  for  next  node  function  in  minimization.  'Hie  cost  of  1  is  determined  by 
d/C(I). 

The  new  don't  care  terms  covered  by  selected  implicant  I  arc  subtracted  from 
d/C(lj.  The  implicant  that  covets  fewer  minteims  has  highci  cost  (is  less  important)  with 
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respect  to  an  implicant  that  covers  many  mintcrms  during  the  minimization  process. 
So  each  don't  care  term  introduced  to  the  remaining  function  affects  A/C(I)  and  in¬ 
creases  its  cost.  Therefore,  MC{ I)  should  be  updated  during  the  minimization  process 


III.  CONSTRAINED  IMPLICANT  SET  HEURISTIC 


A.  INTRODUCTION 

This  heuristic  discussed  in  this  section  is  named  after  the  constrained  implicant  set 
heuristic  (CISH)  because  the  original  concept  has  been  developed  under  the  same  name 
in  multiple- valued  logic  minimization  [Ref.  8].  The  idea  is  to  extend  the  constrained 
implicant  set  concept  in  multiple-valued  logic  to  binary  minimization. 

The  CISH  has  two  computational  phases:  1)  constrain  and  select  a  minterm  ac¬ 
cording  to  its  implicant  cover  size  and  2)  constrain  and  select  an  implicant  with  respect 
to  its  minterm  coverage.  The  selection  of  the  minterm  differs  from  other  exact  minimi¬ 
zation  algorithms  due  to  the  difference  in  its  rules  of  decision. 

The  selection  of  the  implicant  which  covers  a  selected  minterm  depends  on  its  cost 
to  the  minimization.'  CISH  chooses  an  implicant  and  investigates  its  effect  on  the  future 
selection  of  minterms  and  implicants. 

B.  MINIMIZATION  ALGORITHM  OF  CISH 

The  general  steps  of  CISH  are  described  (see  Appendix  G  for  the  C  program  listing) 
in  this  Section.  In  the  algorithm  below,/ denotes  the  function  to  be  minimized. 

All  information  about  minterms,  the  implicants,  and  necessary  sets  are  initialized 
once  and  updated  in  the  recursive  part  of  the  heuristic, 

1.  Initialization 

•  Form  the  uncovered  minterm  set  ( UA1S )  from  all  minterms  of / 

•  Form  the  don't  care  set  (DCS)  from  all  don't  care  minterms  of / (Initially,  this  is 
empty), 

•  Form  the  implicant  set  (IS)  from  all  implicants  of /. 

•  Tally  the  .WC(1)  for  each  implicant  1  in  IS. 

•  Find  the  lCS(o.)  of  each  minteim  in  VMS. 
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2.  Recursive  Algorithm 

Apply  the  following  steps  recursively  to  a  function /until  the  function  consists 
of  only  don't  care  or  0  terms. 

•  Select  an  uncovered  minterm  a  with  the  lowest  7CS(a)  from  UMS  .  If  more  than 
one  such  minterm  exists,  the  one  with  smaller  binary  representation  is  selected. 

•  Construct  the  constrained  implicant  set  C7S(a)  that  includes  all  implica*.ts  that 
•  cover  minterm  a. 

•  If  7CS(a)  #  2,  select  the  implicant  I, (a)  with  the  lowest  MC(I). 

•  If  ICS{«)  «  2  and  A/C(l)  *=  2  for  each  of  the  two  implicants  in  CIS(a),  apply  the 
extended  search  technique  (LSI)  to  select  l,(a), 

•  Put  I, (a)  into  the  minimized  sum-of-products  set  (MSP). 

•  Find  the  minterms  in  VMS  that  are  covered  by  I, (a). 

•  Remove  these  minterms  from  UMS,  and  place  them  into  DCS. 

•  For  each  implicant  in  IS  that  covers  at  least  one  new  don't  care  term,  subtract  the 
number  of  new  don't  care  terms  covered  by  1(a)  from  MC( I).  If  the  A/C(I)  -  0, 
remove  1(a)  from  IS. 

3.  Extended  Search  Technique  (EST) 

Appl>  this  search  technique,  if  there  is  a  minterm  with  7CS(a)  =  2  in  the  recur¬ 
sive  part  of  the  heuristic. 

•  Compute  7CS(/i)  of  each  direct  neighbor  /i  of  a. 

•  Choose  a  /?  with  ICS(p)  =  7CS(a)  (  =  2),  if  such  a  /Sexists.  Select  I, (a)  that  covers  (t 
and  a. 

•  Otherwise,  select  a  neighbor  /?  of  a  with  the  smaller  ICS(p).  Select  I, (a)  that  covers 
/?  and  a. 

C.  EXPLANATION  OF  CISII 

In  this  section,  the  algorithm  described  in  Section  II1.B  is  explained  by  using  ex¬ 
amples. 
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Initialization  Part: 

All  information  about  minterms  and  implicants  are  found  and  prepared  for  the  re¬ 
cursive  phase.  Mainly,  three  basic  working  sets  ( UMS ,  DCS,  IS)  and  two  basic  meas¬ 
ures  ( ICS(ct ),  A/C(I))  of  C1SH  have  been  formed. 

Firstly,  the  UMS  is  formed.  This  set  includes  the  minterms  that  belong  to  A0N  of 
function/  DCS  is  assumed  to  be  empty  at  the  initialization  part.  DCS  and  UMS  are 
complimentary  sets.  When  UMS  decreases,  DCS  increases.  The  IS  initially  includes  all 
implicants  of  the  function.  IS  gets  smaller  as  a  result  of  some  implicants  being  removed 
during  the  minimization  process. 

Two  basic  measures  of  C1S11  are  formed  in  initialization  part.  These  are  A/C(I)  and 
ICS( a).  All  minterm's  7C5(a)  are  computed.  The  CIS1I  computes  the  /CS(a)  the 
column-row  older.  For  example,  the  minterm  a.  -  1  (i.c.,  binary  representation  0001)  is 
evaluated  earlier  than  minterm  a**  2  (i.e.,  binary  representation  0010). 

A/C(I)  is  calculated  by  counting  the  number  of  the  minterms  covered  by  I.  For  example, 
if  1  covers  4  minterms,  its  A/C(l)  is  equal  to  4.  But.  the  A/C(l)  of  each  1  will  change 
during  the  minimization  process  as  mentioned  in  Section  I  I.C. 

Recurs  he  Part: 

Recursive  part  is  the  second  computational  phase  of  C1S11.  The  purpose  of  the  re¬ 
cursive  part  is  to  select  the  minterm  with  the  lowest  ICS( a),  then  constrain  and  select  the 
implicant  that  has  the  lowest  cost  to  the  minimization.  All  of  the  minimization  has  been 
done  in  this  part  of  heuristic.  All  the  steps  shown  in  Section  III.B  arc  explained  by 
Fxample  6  and  7  in  this  section. 

It  is  important  to  select  the  fust  minterm  intelligently  during  the  minimization 
process  in  C1SH.  The  importance  of  the  nrinterms  is  determined  by  ICS{a),.  As  men¬ 
tioned  in  Chapter  II.  ICS(c/)  is  a  measuic  of  how  many  possible  combinations  a  minteim 
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has  with  neighboring  terms.  A  lower  lCS(a)  means  that  fewer  combinations  exist.  The 
minterm  with  lower  1CS(&.)  gains  importance  with  respect  to  the  minterms  having  higher 
/CS(ot)  in  CISH,  because  the  minterm  with  lower  ICS( a.)  tends  to  be  isolated  and  the 
minterm  with  a  higher  ICS(a)  tends  to  be  in  the  middle  of  the  cluster  of  1  's., 

If  a  function  has  £  minterms,  the  search  space  of  the  function  has  t,  main  branches 
(sec  Figure  11  in  Section  D).  C1SH  selects  the  minterm  with  the  lowest  /CS(a)  in 
VMS,  then  selects  the  main  branch  that  includes  the  selected  minterm.  All  of  the  min¬ 
imization  process  in  C1SH  is  done  in  this  main  branch  and  its  subbranches.  CISH  does 
not  make  a  search  to  find  the  minimal  expression  of  the  function  for  any  one  of  the  re¬ 
maining  £-1  main  branches.  One  of  the  properties  of  CISII  is  to  find  a  near  minimal 
solution  by  searching  only  one  main  branch  of  over-all  search  space.  That's  why  the 
selection  of  the  minterm  is  very  important.  The  selection  of  minterm  is  done  at  every 
node  function  generated  from  root  function  by  heuristic  rules.  The  selection  of  minterm 
with  the  lowest  ICS{ a)  reduces  the  search  space  significantly  (unlike  the  exact  minimi¬ 
zation  algorithms  such  as  Quine- McCluskev), 

Example  5: 

'1  o  illustrate  the  effects  of  the  selection  of  minterms  in  minimization,  consider  the 
following  four  variable  function  with  eight  minterms. 

Let 


J[w,  x,y,  2)  =  V  in(  (j,  1, 2,  3,  4,  7,  9,  10  )  ,■ 

Two  possible  minimal  solutions  for  the  given  function  are  shown  in  Figure  7  on  page 
23. 
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If  the  minimization  begins  at  any  one  of  the  mintcims  in  the  center  of  a  cluster  (i.e., 
a  =  0,  1,  2,  3),  the  prime  linplicant  that  corresponds  to  ir.v  will  be  included  in  the  solution 
sct.;  This  prime  implicant  is  indicated  by  a  dashed  line  in  Figure  7. a.  Four  additional 
implicants  aic  necessary  to  co\er  the  remaining  minterms  with  ICS{ c)=  1.  The  sum-of- 
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products  expression  of /  includes  five  product  terms.  The  implicant  wx  is  redundant  in 
the  minimization  of f 

If  the  minimization  begins  with  any  one  of  the  minterms  such  that  lCS(a)=  1  (i.e., 
a»»  4,  7,  9,  10),  the  minimal  solution  will  have  four  prime  implicants  as  seen  in 
Figure  7.b.  The  sum-of-products  expression  of / for  two  solutions  are  written  below: 

J[w,  x,y,  z)  -  wyz  +  xyz  +  wyz  +  xyz  +  wx 
J[w ,  jc,  j/,  z)  -  wyz  +  xyz  +  wyz  +  xyz 

The  CISH  constructs  the  C7S(«)  after  selection  of  minterm  from  the  implicants  in 
IS.  A  group  of  implicants  in  IS  is  constrained  by  using  ICS( a).  All  implicants  in  C/S(a) 
cover  the  minterm, a. 

The  selection  of  implicants  is  equivalent  to  breaking  the  coupling  between  that 
implicant  and  its  neighbors.  The  candidate  implicant  should  have  the  lowest  cost  to 
minimization  in  CISH.  [Ref.  9] 

The  coupling  strength  is  introduced  as  the  minimization  cost.  The  A/C(I)  is  used  as 
the  cost  of  implicant  to  minimization.  The  implicant  with  the  lowest  cost  has  the  highest 
A/C(I).  The  implicant  with  the  highest  A/C(l)  covers  the  largest  area  in  the  function. 
The  largest  area  can  contain  don't  care  terms  as  well  as  /  terms. 

The  chosen  implicant  in  C/S(a)  is  to  be  a  prime  implicant  due  to  the  prime  implicant 
theorem.  The  prime  implicant  theorem  states  that  a  minimal  sum-of-product  must  al¬ 
ways  consist  of  a  sum  of  prime  implicants.  [Ref.  10:  p.  206] 

At  least  one  of  the  implicants  in  CIS(a)  should  be  a  prime  implicant  and  it  contains 
the  fewest  literals  among  other  implicants  that  cover  minterm  a.-  Naturally,  the  prime 
implicants  always  have  higher  A/C(I)  than  non-prime  implicants  in  C/S(a),  The  strategy 
of  CISH  is  to  always  select  an  implicant  that  covers  minterm  a  with  the  lowest  cost.; 
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This  rule  provides  for  that  selected  implicant  being  one  of  the  prime  implicants  having 
the  highest  A/C(I)  in  C/S(a). 

Extended  Search  Technique: 

This  technique  is  applied  only  for  the  special  case  when  the  /CS(a)  of  a  minterm  a 
equals  two.  It  indicates  that  minterm  a  is  covered  by  only  two  implicants  and  at  least 
one  of  them  should  be  chosen  as  a  part  of  the  solution.  If  the  A/C(I)  of  these  implicants 
is  not  equal  to  two,  CISII  selects  the  implicant  with  the  highest  A/C(l).  On  the  other 
hand,  if  the  A/C(I)  of  each  implicant  in  CIS( a.)  equals  two,  it  indicates  that  the  minterm 
a  is  covered  by  two  implicants,  that  none  of  them  covers  don't  care  terms,  and  their  cost 
to  the  minimization  process  is  equal.  In  this  case,  the  best  and  the  most  efficient 
implicant  should  be  chosen  in  ClS(x),  The  purpose  of  the  EST  is  to  minimize  the  neg¬ 
ative  impact  for  future  minterm  selection  as  well  as  implicant  selection  by  choosing  the 
most  efficient  implicant  in  C7S(a). 

In  EST,  CISII  finds  the  two  direct  neighbors  /?  of  minterm  a.  It  checks  the  ICS(P) 
of/?.  If  a  P  exists  such  that  ICS{P)  =  ICS(<x)  ( =  2),  it  selects  the  implicant  in  C/S(a)  that 
covers  both  /?  and  a.  If  both  direct  neighbors  have  the  same  1CS(P)  such  that 
ICS(p)=  ICSia),  CISII  selects  the  neighbor  with  smaller  binary  representation,  then 
chooses  the  implicant  that  covers  both  P  and  a.  If  none  of  its  direct  neighbors  have  the 
same  JCS(p)  as  ICS{ a).  CISII  then  selects  the  P  with  smaller  ICS(P). 
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Example  6: 


To  illustrate  the  application  of  EST  in  CISH,  consider  the  following  a  four  variable 
function  with  eight  minterms.  Options  for  the  minimization  of / will  either  1)  not  apply 
EST  (Figure  8.a)  or  2)  apply  EST  (Figure  8.b)  as  shown  on  page  27. 

Let 


J[w,x,y,z)**  ^  m{  2, 3, 5, 6, 7,  9, 11, 13). 

There  arc  four  minterms  (i.e.,  a  —  5,  9,  11,  13)  to  be  selected  due  to  their  ICS{u)  (i.c., 
2).  The  CISH  selects  the  minterm  a  =  5  (located  at  0101)  according  to  column-row  or¬ 
der,  The  selected  minterm  has  two  implicants  that  cover  it  (01-1,  -101),  that  is, 
CISfOlOl )  has  implicants  I, (a)  and  I2(a)  represented  in  binary  as  01-1  and  -101  respec¬ 
tively. 

Figure  S.a  shows  the  case  when  EST  is  not  applied.  Since  two  implicants  (!,(«), 
J2(oc))  hate  the  same  cost  for  minimization  (i.c.,  MC(I)  =  2),  the  one  of  them  must  be  se¬ 
lected.  If  l,(cO  is  selected  as  the  implicant  to  be  placed  in  the  solution  set,  then 
Figuie  S.a  shows  the  minimal  sum-of-products  expression  of  the  function. 

On  the  other  hand,  if  EST  is  applied,  after  selection  of  the  minterm  a  =  5,  its  direct 
neighbors  arc  found  as  /?,  =  01 1 1  and  /?2=1111.  Since  1CS!0111)~A  and 
ICSnill)-  2,  EST  selects  1111  which  has  lower  value.  In  words,  the  implicant  I2(a) 
(that  covers  the  a  and  /?2)  is  the  best  and  the  most  efficient  implicant  to  minimize  nega¬ 
tive  impact  for  future  selection  of  minterms  and  implicants  at  the  next  node  function. 

It  can  be  seen  that  if  EST  is  not  applied,  there  may  be  a  negative  impact  to  min¬ 
imization  of  function.  The  minimal  solution  obtained  in  this  way  has  four  product 
terms.  EST  provides  a  better  minimal  solution  for  the  same  function  .  The  minimal 
solution  has  only  three  product  terms  as  shown  in  Figure  S.b. 
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Figure  9.  Constrained  linphcants  in  Extended  Search  technique 

The  selected  mintenns  and  implicants  in  binary  representation  is  placed  next  to  these 
lines.  Another  solution  (without  applying  LSI)  is  indicated  by  italic  font  located  on  left 
branch  front  the  root.; 

The  selected  implicant  l,(a)  is  placed  into  MSP.  C1SH  finds  the  mintenns  covered 
by  I, (a)  in  UMS.  It  turns  these  mintenns  to  the  don't  care  terms  and  puts  them  into 
DCS.-  CISI1  updates  the  cost  of  the  implicant  in  IS  that  covers  at  least  one  of  the  new 


don't  care  term.  Updating  of  the  cost  of  1(a)  is  done  by  subtracting  the  number  of  newly 
introduced  don't  care  terms  (covered  by  1(a))  from  A/C(I). 

For  any  implicant  I,  if  A/C(I)  =  0,  then  it  is  removed  from  IS,  The  number  of 
implicants  in  the  IS  gets  smaller  in  each  recursive  step.  The  computation  stops  when 
VMS  is  empty. 

D.  A  WALKTHROUGH  EXAMPLE  OF  CISH 

It  is  instructive  to  examine  the  application  of  the  CISH.  A  four  variable  function 
with  eight  minterms  is  used  as  an  example.  The  input  function  has  been  shown  in 
Figure  10  on  page  30. 

Let 


x,y,  z) «  £  m(  3. 4, 5, 7, 9. 13. 14, 15 ). 

All  information  about  minterms  and  implicants  are  presented  in  Table  1  and  Table  2 
on  page  32.  For  simplicity,  binary  representation  of  cadi  implicant  is  indicated  by  a 
capital  letter  in  Table  1, 

The  sets  of  CISH  constructed  in  the  initialization  part  are  listed  below: 

VMS  —{3,4,  5,  7,  9, 13, 14,  15  } 

IS  =  {  A( 2),  B( 2),  C(2),  D{2),  £(4),  F( 2),  G(2),  1(2),  J( 2) } 

DCS  =  4> 

CIS(a)  =  </> 

There  aie  eight  main  branches  in  the  search  space.  These  branches  are  shown  in 
figure  11  on  page  34,  Cl  SI  1  can  begin  minimization  from  any  one  of  these  four 


branches  (a  =  3,  4,  9,  14),  because  their  /CS(«)  are  the  same  and  the  lowest  in  the  VMS 
(i.e.,  ICS{oi)=  1),  By  using  ICS( a),  CISH  reduces  the  search  space  from  eight  branches 
to  four  branches.. 

It  is  assumed  that  mi'ntcrm  a  =  3  is  selected  to  begin  minimization  among  four 
minterms  with  ICS  (u)  -  2.  CIS(c/)  is  constructed  with  respect  to  o.  =  3.  This  set  consists 
of  only  i'mplicant  I)  in  IS.  Implicant  I)  covers  a=3  and  a  =5.  The  sets  of  C1S11  aic 
updated  with  the  rules  of  the  heuristic.;  For  example,  the  costs  of  implicant  E,  F,  I 
change  in  IS.  The  cost  of  each  implicant  is  indicated  in  parenthesis  next  to  the 
implicant.  The  new  sets  of  CIS1I  and  costs  of  the  implicants  are  shown  below  and  cor¬ 
responding  to  the  leftmost  node  at  level  1  in  Figure  11  on  page  34, 

VMS  =  {4,5,  9,  13,14,15} 
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IS  =  {  A( 2),  C(2),  D( 2),  E( 3),  F(t),  (7(2),  7(1),  7(2) } 
DCS  *={3,7} 

MPS  =  {  B } 


Table  1.  INFORMATION  ABOUT  IMPLICANTS 

Implicant  (1(a)) 

Binary  Represen¬ 
tation 

A/C(I) 

M  interm  Covered 
by  1(a) 

A 

010- 

2 

4,5 

B 

0-1 1 

2 

3,7 

Tahle  2.  INFORMATION  ABOUT  MINTERMS 


UMS  includes  six  minterms  and  CIS11  chooses  three  of  them  to  start  the  minimiza¬ 
tion,  since,  their  ICS{a)  is  smaller  than  others.  These  are  a  =  4,  9,  14  and  their  ICS{«) 
are  equal.  It  is  assumed  that  a- 4  is  selected.  C/S(a)  contains  only  implicant  A.  A 
covers  a  =  4  and  5.  The  new  sets  of  CISH  becomes: 

VMS  =  {9,13, 14,15} 

IS  -  {  C(2),  D{2),  E{2),  G(2),/(1),J(1)} 

DCS  m  {  3, 4,  5,  7  } 

MPS  -  { A,  D) 

Implicant  F  was  removed  from  IS  because  it  would  cause  the  highest  cost  for  all  future 
minimization  processes  (i.e.,  A/C(I)  =  0 ). 

There  are  only  two  minterms  to  begin  the  next  selection  and  constraining  implicant. 
These  are  a  =  9.  14  and  they  arc  equal  in  ICS(a).  So  we  can  arbitrarily  select  the  one 
with  smaller  binary  representation,  i.e.,  minterm  a  =  9  is  selected.  The  C/S(a)  includes 
only  implicant  C. 

After  selection  of  C,  the  sets  are  updated  as  shown  below:. 

VMS  =  {  14, 15  } 

IS  =  {  D{2),  £(!)..  G(l),  7(1) } 

DCS  =  {3,4,  5,  7,  9,  13  } 

MPS  -  {A.  D,C } 
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Now,  only  one  minterm  is  in  UMS  (  or=  14).  CIS{ a.)  consists  of  only  implicant  D. 


D  covers  a=  14  and  a  -  15.  The  sets  of  the  CISI1  become  : 
UMS  =  <p 


DCS  ={3,4,  5,  7,9,13, 14,15} 

MPS  *  { A,  B,C,D] 

The  minimal  sum-of-products  expression  of/is  the  UMS  or: 

J{  iv,  x,y,  z  )  *  010  — h0  —  114-1  —  01-4-  111  — 

*=»  vvxv  +  wyz  +  wyz  +  wxy 

All  branches  and  constrained  implicunts  arc  shown  in  Figure  1 1  on  page  34.  The 
constrained  paths  are  shown  with  bold  lines,  located  on  the  left  branch  from  root.  An¬ 
other  search  path  is  shown  by  selection  a2'  5  at  the  very  beginning  of  the  minimization. 
This  path,  located  on  the  right  branch  from  the  root,  gives  another  possible  solution, 
without  applying  CISII, 
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IV.  COMPARATIVE  RESULTS 


In  this  thesis,  all  test  results  were  obtained  by  running  sample  functions  on  the  VAX 
11/785  and  ISI  workstations.  Different  number  of  sample  functions  (4100)  were  ran¬ 
domly  generated.  The  input  functions  are  generated  for  different  variables  with  different 
number  of  minterms,  (i.e.,  9  variable  function  with  475  minterms  or  7  variable  function 
with  120  minterms).  Each  algorithm  was  applied  to  these  sample  functions,  then  the 
average  number  of  product  terms,  average  computation  time,  and  average  memory  usage 
are  recorded.  The  computation  time  for  9  variable  functions  and  larger  is  very  large. 
This  explains  why  we  did  not  simulate  more  than  9  variables.  This  thesis  investigated 
three  algorithms:  1)  CISH,  2)  Mill  (see  Appendix  B),  3)  Espresso  II  (see  Appendix  C). 

A.  PERFORMANCE  COMPARISON 

The  performance  measures  arc  recorded  and  compared.  These  are  1)  the  average 
number  of  product  terms,  2)  the  average  computation  time,  and  3)  the  average  memory 
usage. 

The  average  number  of  product  terms  will  show  us  the  advantage  of  the  each  algo¬ 
rithm.  T  or  each  set  of  sample  input  functions,  the  average  number  of  product  terms  (see 
Appendix  D)  is  computed.  From  these  data  points,  a  curve  is  plotted  to  indicate  the 
as  erage  number  of  product  terms  as  a  function  of  the  number  of  the  variables.  The  plot 
is  shown  m  Figure  12  on  page  36. 

In  this  exponential-growing  shaped  figures,  it  is  observed  that:' 

•  The  differences  in  the  number  of  average  products  terms  among  these  algorithms 
are  not  significant.  Less  than  1%  of  testing  functions  are  different. 

•  When  the  number  of  variables  gets  larger,  the  curse  grows  up  exponentially.  This 
is  chiefh  due  to  the  fact  that  both  the  number  of  minterms  and  nnpheants  increase 
exponentially  and  then  make  the  computation  time  longc.. 


Figure  12.  A\erage  Number  of  Product  Terms 

From  Figure  12,  the  CIS! I  performs  as  well  as  the  other  two  algorithms  (Mill, 
Espresso  II). 

B.  TIMING  COMPARISON 

Both  VAX  11/785  and  IS1  woikstations  at  NFS  can  measure  a  program's  compu¬ 
tation  time  in  UNIX  environment.  In  this  section,  timing  comparison  counts  the  aver¬ 
age  computation  time  for  each  group  of  input  functions. 
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Figure  13.  Average  Computation  lime 

The  decision  rules  for  selecting  the  minterm  n  and  the  implieant  1(a)  in  each  algo¬ 
rithm  are  different.  Generally,  an  algorithm  that  has  complex  rules  to  select  a  and  1(a) 
takes  longer  computation  time.  The  C1SII  uses  more  complex  decision  rules  than  Mill 
and  Espiesso  11.  It  may  appear  that  CISil  needs  mote  computation  time  than  other 
algorithms.  However,  the  computation  time  of  the  C1SI1  is  shorter  than  Mill,  although 
it  is  not  as  last  as  Espresso  II.  'I he  graphical  lcsult  is  shown  in  Eigutc  13.  Espresso  11 
outperforms  the  other  two  hcuiistics.  This  might  be  due  to  smaller  constant  in  compu¬ 
tation  complexity.  The  numciical  results  ate  shown  in  Appendix  E. 


C.  MEMORY  COMPARISON 

The  average  memory  usage  is  measured  for  each  algorithm  (see  Appendix  F).  From 
these  data  points,  a  curve  is  plotted  to  indicate  the  average  memory  usage  as  a  function 
of  the  number  of  the  variables..  The  plot  is  shown  in  Figure  Id  on  page  38. 

It  is  expected  that  a  heuristic  or  algorithm  based  on  tabular  method  needs  more 
memory  to  store  the  information  for  keeping  track  of  each  term  during  the  minimization 
process.  The  Mill  is  based  on  a  tabular  method,  thus  it  needs  a  lot  of  memory  space. 
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The  CISH  selects  and  constrains  the  minterm  a  and  implicant  1(a).  As  a  natural  con¬ 
sequence,  CISH  uses  less  memory  than  MIH  due  to  heuristic  strategy. 

A  sophisticated  dynamic  memory  scheme  was  used  in  Espresso  II,  so  that  it  de¬ 
mands  less  memory  than  the  other  two  heuristics  (MIH,  CISII)., 

To  summarize: 

•  There  has  been  significant  difference  among  the  memory  usage  of  each  algorithm. 

•  Direct  covering  and  tabular  method  used  in  CISH  and  MIH  needs  more  room  to 
store  the  information  about  terms  than  the  decomposition  technique  used  in 
Espresso  II. 

•  If  we  use  dynamic  memory  allocation  in  CISH,  the  memory  requirement  will  be 
reduced. 


V.  DISCUSSIONS  AND  CONCLUSIONS 


A.  DISCUSSIONS 

There  is  no  significant  difference  among  these  algorithms  with  respect  to  number 
of  product  terms.  The  small  difference  in  the  average  number  of  product  terms  origi¬ 
nated  from  the  application  of  EST  to  CISH.  Recall  that  in  EST,  we  compute  /CS(ot) 
of  a  given  minterm,  then  constrain  and  select  the  implicant  I, (a)  with  respect  to  the  re¬ 
lationship  among  direct  neighbors  of  minterm  a  (i.e,,  one  step  look-ahead). 

Naturally,  a  further  look-ahead  while  selecting  the  I, (a)  may  provide  a  better  se¬ 
lection  of  implicants.  The  exponential  growth  of  the  number  of  all  possible  implicants 
restricts  the  practical  use  of  k  look-aheads  for  k  +  1.  The  application  of  EST  provides 
better  solutions  for  CISH  over  MIH.  On  the  other  hand,  the  best  solutions  are  provided 
by  Espresso  II. 

From  the  computation  time  results,  the  CISH  runs  faster  than  MIH  in  all  testing 
conditions.  This  time  efficiency  is  a  result  of  the  decision  rules  employed  in  CISH  (that 
takes  advantages  of  the  properties  of  the  constrained  implicant  set  concept).  On  the 
other  hand,  Espresso  II  runs  faster  than  the  other  two  heuristics  as  a  results  of  the  be¬ 
havior  of  the  algorithm  based  upon  the  unate  paradigm  (see  Appendix  C).  Typical 
PLA's  produce  shallow  recursion  trees  terminating  quickly  at  unate  leaves  in  search 
space.  This  benefit  is  used  in  Espresso  II. 

The  memory  comparison  shows  that  MIH  and  CISH  need  more  room  to  store  the 
information  about  the  minterms  and  the  implicants  to  operate  efficiently..  On  the  other 
hand,  Espresso  II  uses  less  memory  by  dealing  primarly  with  matrix  representation  of 
the  function  and  the  minimization.  It  is  believed  that  use  of  more  memory  makes 
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heuristics  slower,  since  memory  intensive  programs  tend  to  run  slower  than  cpu  intensive 
programs. 

B.  CONCLUSIONS 

One  of  the  primary  goals  of  this  research  is  to  examine  whether  the  constrained 
implicant  set  concept  in  multiple-valued  logic  can  be  efficiently  used  in  binary  minimi¬ 
zation.  It  is  seen  that  the  binary  minimization  of  a  given  function  can  be  done  efficiently 
by  using  the  constrained  implicant  set  concept. 

In  the  development  of  the  C1SII  and  from  comparative  results,  we  have  the  follow¬ 
ing  observations: 

•  The  constrained  implicant  set  concept  reduces  the  search  space  significantly  in  the 
binary  minimization. 

•  C1SH  does  not  lose  run  time  efficency  unlike  Mill  because  the  heuristic  finds  the 
solution  and  stops  earlier  than  Mill. 

•  Direct  covering  fits  well  with  the  constrained  implicants  set  concept.  On  the  other 
hand,  it  does  not  provide  efficiency  on  the  computation  time  as  much  as  the  de¬ 
composition  technique  docs,  such  as  Espresso  II. 

•  In  the  cyclic  case,  applying  EST  provides  better  results  than  Mill.  On  the  other 
hand,  using  onh  one  step  look-ahead  in  ESI'  may  loose  the  optimality.  However, 
it  is  not  practical  to  have  mote  than  one  look-ahead,  that  requires  longer  compu¬ 
tation  time  and  larger  memory  space.. 

•  By  constraining  implicants  with  the  cost  factor,  it  is  possible  to  find  near  and  good 
minimal  sum-of-products  expressions. 

•  It  is  possible  to  get  a  near  and  good  minimal  solution  by  only  searching  the  main 
branches  of  the  search  space. 

•  The  memory  usage  can  be  decreased  by  using  the  dynamic  memory  allocation  like 
Espresso  II.  Originally  OSH  uses  the  static  memory  allocation. 


Generally,  Espresso  II  is  more  efficient  than  CISH.  This  is  not  surprising  since  our 
heuristic  is  an  initial  unoptimized  prototype.  Espresso  II  on  the  other  hand  is  the 
product  of  considerable  effort  b\  a  very  large  team  over  a  long  period  which  has  con¬ 
centrated  on  producing  a  production  quality  software  package. 


Directions  for  Future  Research 

It  is  our  hope  that  the  ideas  and  heuristic  described  in  this  thesis  represent  a  plateau 
for  the  development  of  two-level  binary  minimization.  The  minimization  problem  is  in 
principle  difficult,  and  future  developments  will  exploit  fundamentally  new  ideas.  Here, 

we  briefly  describe  directions  in  which  future  research  might  be  pursued, 

•  In  C1SH,  function  is  represented  in  the  form  of  the  truth  table.  All  the  minimiza¬ 
tion  process  is  done  by  using  the  data  structure.  On  the  other  hand,  Espresso  II 
uses  the  matrix  representation  of  function  and  the  minimization.  This  provides 
increased  speed  of  execution.  This  matrix  representation  can  be  applied  CISH. 

•  Direct  covering  technique  is  applied  to  CISH,  It  is  obtained  that  decomposition 
technique  can  provide  faster  and  more  accurate  solutions  in  minimization  like 
Espresso  II.  The  constrained  implicants  set  concept  can  be  improved  by  using  the 
decomposition  technique  instead  of  direct  covering. 

•  The  speed  of  execution  can  be  increased  by  applying  the  concept  of  unate  paradigm 
to  CISH.  Actually,  Cl  SI  I  reduces  the  search  space  significantly.  On  the  other 
hand,  we  believe  that  using  the  concept  of  unate  paradigm  reduces  the  search  space 
(that  already  reduced  by  applying  the  constrained  implicant  set  concept). 

•  The  function  applied  to  CISH  includes  only  the  1  and  0  terms.  The  CISH  can  be 
improved  by  using  the  don't  care  terms  at  the  beginning  of  the  minimization.  We 
hope  that  using  the  don't  care  terms  at  the  beginning  of  the  minimization  with  1 
terms  helps  reduce  the  computation  time  and  use  less  memory  space. 

•  CISH  is  developed  using  only  the  single  output  case.  The  heuristic  can  be  modified 
to  be  used  in  the  minimization  of  the  multi-output  case.  We  believe  that  CISH 
provides  more  efficient  results  than  Mill  in  the  multi-output  cases. 
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APPENDIX  A.  ABSOLUTE  MINIMIZATION  ALGORITHM 


This  algorithm  is  taken  from  [Ref.  1]. 


***  ft***  #»*$***#  ft*#*'******************#*’!'**** 


algorithm  absolutejninimization; 
f  ♦-  input_function; 

cur_best_soln_sct  ♦-  best  solution  from  the  Fompcr  and  Armstrong, 
Besslich  and  Ducck  and  Miller  heuristics; 
cur_bcst_soln_sizc  «-  number  of  implicants  in  best  solution; 
cur_partial_soln_set  ♦-  «/>; 
cur_partial_soln_sizc  0  ; 

minimized  0; 

stop 


procedure  minimize(f); 

<f>  some  essential  implicant  set  of  f ; 

while  (( there  exists  another  implicant  in  <p)  and  (  cur_partial_soln_si/e  +  1  < 

cur_best_soln_si/.e  ) )  do 

begin 

I  <-  the  next  implicant  in  <!> 
cur_partial_so!n_set  <-  cur_paitial_soln_set  (J{I) 


4.' 


cur_partia!_soln_$ize  ♦-  cur_partial_soln_size  +  1; 

if(  for  all  assignments  x  of  values  to  X,  =  0  or  f(x)  =  r ) 

then 

begin 

cur_bcst_soln_set  «-  cur_partial_soln_set; 
cur_best_soln_set  cur_partial_soln_size 

end 

[[backtrack  ]] 

else  if  (  cur_partial_soln_size  +  1  <  cur_best_soln_size  )  then 
minimize(  f  *  I); 

cur_partial_soln_size  <-  cur_partial_soln_size  -  1; 
cur_partial_sohi_set  «-  cur_partial_soln_sct  -  {1} 
end 
return 


.>11(1*  ijl*****!!.*!!!***********!)'******!!'***'****’!"****** +****  +  *  +  *  e*************** 


;*  The  subtraction  of  an  iinplicant  1  from  a  function  f,  as  described  by  f «-  f  - 1, 
takes  into  account  the  value  of  the  input_function.  Let  x  be  assignment  of 
values  to  variables  X,  Then,  f  <-  f  - 1  means 
for(all  assignments  x  of  values  to  X)  do 
begin 

if((  f(x)  =  r  )  or  ( input_  function(x)  =  r )  or  (l[x)  <  I  (x)  and  input_function(x) 

=  r-1)) 


then  f(x) «-  r 

else  f(x)  *-  f(x)  -  I  (x) 
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end 


I 

*********  +  ^j>*$**^  ****************************  *****  +  ****************** 
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APPENDIX  B.  MAXIMUM  IMPLICANT  HEURISTIC 


The  Maximum  Implicant  Heuristis  (MIH)  is  developed  by  modifying  the 
Quine-McCluskey  Algorithm  (QM).  MIH  is  based  on  obtaining  the  switching  function 
as  a  near  minimal  sum-of-products.  MIH  differs  from  QM  with  respect  to  the  two  basic 
approaches.  These  are  stated  below: 

•  Incorporation  of  A0N,  AorF,  ADC  sets  of  function  at  the  beginning  of  the  minimiza¬ 
tion 

•  Approaching  to  solve  the  cyclic  case  for  a  given  function 

The  differences  between  two  algorithms  are  explained  by  showing  the  main  steps  of 
QM,  QM  is  an  exact  minimization  algorithm.  It  consists  of  two  main  parts  shown  be¬ 
low  [Ref.  2::  P.  58]: 

•  Generation  of  the  prime  implicants 

•  Extraction  of  a  minimum  prime  cover 

Modifications  are  done  in  these  main  parts  of  QM. 

Modification  1: 

The  first  modification  to  develop  MIH  is  done  at  the  first  part  of  QM,  The  Mill 
does  not  include  the  don't  care  terms  while  generating  all  prime  implicants.  It  uses  only 
the  minterms  that  belong  to  Aav  and  AofF  of  given  function  at  the  beginning.  On  the 
other  hand,  QM  can  consider  Aos,  Aorr,  ADC  of  function.  Use  of  the  don't  care  terms  in 
the  minterm  list  makes  generation  of  all  the  prime  implicants  more  complex,  The  more 
computation  time  and  memory  space  is  necessary'  to  find  prime  implicants.  Thus,  a 
modification  has  been  done  to  save  the  computation  time  and  memory  space.  This 
modification  is  reflected  to  de\elopment  of  MIH  as  not  including  ADC-  of  a  gi\en  function 
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at  the  beginning  of  minimization  and  dealing  with  only  the  function  with  ACN  and  A0FF 
in  its  minterm  list. 

Modification  2: 

QM  algorithm  sets  up  a  covering  table  that  shows  all  of  the  prime  implicants  and 
mintcrms  covered  by  prime  implicants.  The  essential  prime  implicants  are  removed  from 
this  covering  table,  then  a  reduced  covering  table  is  formed.  Secondary  essential  prime 
implicants  are  found  and  removed  from  last  covering  table.  A  new  covering  table  is  set 
up  by  remaining  minterms  and  prime  implicants  [Refs.  2:  pp.  59-64,  11:  pp.  146-156]. 

In  fact,  it  occurs  in  last  reduced  covering  table,  there  is  more  than  one  possible  cover 
for  given  function.  It  is  said  that  function  has  cyclic  case.  This  cyclic  case  is  solved  by 
apllying  the  Petrick  Algorithm  in  QM.  Petrie  Algorithm  can  be  stated  as  producing  all 
possible  covers  for  the  function  in  covering  table  and  selection  of  the  one  of  the  covers 
requiring  the  smallest  number  of  prime  implicants  and  literals  [Ref.  2:  p.  64]. 

Producing  all  possible  covers  for  the  function  needs  more  computation  time  and 
memory  requirement.  All  possible  co\ers  in  last  reduced  covering  table  must  be  found 
in  QM.  The  reason  is:  QM  quarantecs  exact  minimal  solution  by  searching  all  possible 
covers  and  selecting  one  having  the  fewest  number  of  product  terms. 

The  second  modification  is  made  at  this  part  of  QM«  The  reason  is  to  avoid 
spending  a  lot  of  computation  time  to  find  all  possible  covers  and  using  more  memory 
place. 

The  modification  can  be  stated  as:'  1)  find  the  maximum  implicant  that  covers  the 
largest  area  in  the  last  reduced  covering  table  and  2)  declare  it  being  in  the  solution  set. 
The  minimization  process  is  continued  after  remosing  the  maximum  prime  implicant 
from  this  co\ering  table,  then  establishing  a  new  reducing  table.  If  function  still  has  a 
cyclic  case,  the  same  process  is  applied  until  no  cyclic  case  exists. 
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The  example  about  applying  Petrick  Algorithm  in  QM  and  breaking  the  cyclic  case 


in  MIH  can  be  found  in  (Ref.  2:  pp.  64-65]. 


Minimization  Algorithm  of  iMlH: 

In  summary,  the  MIH  (for  finding  a  near  minimal  sum-of-products  expression  to  a 

given  function)  follows  the  step  given  below. 

•  Find  the  set  of  all  prime  implicants  of  the  function  by  using  minterms  that  belong 
to  A0N  and  A0FF  of  the  given  function. 

•  Construct  a  covering  table  from  all  generated  prime  implicants. 

•  Identify  all  of  the  essential  prime  implicants  and  form  a  reduced  covering  table. 

•  Identify  secondary  essential  prime  implicants  and  reduce  the  covering  table  again. 

•  If  there  is  a  cyclic  case  in  reduced  covering  table,  break  the  cyclic  case  by  selecting 
the  maximum  prime  implicant  to  find  a  minimal  cover  for  remaining  minterms. 
Apply  this  process  until  no  cyclic  case  occurs  in  reduced  covering  table. 


Observations  on  MIH: 

The  most  important  observations  about  MIH  are  stated  below: 

•  Most  of  the  steps  of  both  algorithms  (MIH,  QM)  are  identical.  If  a  function  does 
not  have  any  cyclic  case  during  its  minimization  process,  Mill  and  QM  can  be 
considered  as  the  same  algorithm. 

•  The  MIH  provides  a  near  minimal  sum-of-products  expression  for  given  function 
as  a  results  of  not  finding  all  possible  covers  to  solve  cyclic  case  in  reduced  cover¬ 
ing  table.; 
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APPENDIX  C.  ESPRESSO  II 


Espresso  II  is  a  set  of  algorithms  for  logic  minimization  which  basically  follows  the 
sequence  of  top-level  transformation  of  iterated  expansion-reduction  pionereed  by  MINI 
l Ref.  5:  p.13]. 

The  fundamental  definitions  used  in  explanation  of  Espresso  11  are  summarized  be¬ 
low: 

Definition : 

A  logic  function  /  is  monotone  increasing  (monotone  decreasing)  in  a  variable  x,  if 
changing  x )  from  0  to  1  causes  all  the  outputs  of /  that  change,  to  increase  also  from  0 
to  1  (from  1  to  0).  A  function  that  is  either  monotone  increasing  or  monotone  de¬ 
creasing  in  x,  is  said  to  be  monotone  or  unate  in  .v,. 

Definition  : 

A  function  is  said  to  be  unate  function,  if  it  is  unate  in  all  its  variables.  Example: 

/=  *i-Vj  4 


Theorem : 

'l  he  Shannon  expansion  theorem  states  that  a  function  can  be  expanded  about  any 
chosen  variable  that  is,  to  produce  an  equivalent  expression  for  the  function  in  which  the 
chosen  variable  appears  once  in  uncomplemented  form  and  once  in  complemented  form. 
The  statement  of  the  theorem  is: 


Jix^.x^tX,,  ».,xn)  =  xtf[x].xi,  1,  .,,,xn)  +  VUv-Vj,  €..0, 
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Algorithms  of  Espresso  II  based  on  a  single  basic  strategy:  a  recursive  divide  and 
conquer.  It  basically  uses  decomposition  technique.  Decomposition  is  based  upon 
Shannon  Expansion.  The  Shannon  Expeansion  uses  the  cofactors  of  a  logic  function. 
Since  Espresso  II  uses  the  benefit  of  unate  functions  in  Shannon  expansion. 

The  algorithms  of  Espresso  II  form  a  logic  minimization  tools  which  actieves  both 
robust  performance  and  quality  results.  Iterative  improvement  produces  well-minimized 
cover  with  high  confidence,  while  unite  paradigm  together  with  special-case  handling 
ensures  reasonably  efficient  execution  for  a  broad  range  of  incoming  problem  [Ref.  5: 

pp.  161-162]. 

Speed  of  execution  in  Espresso  II  is  based  on  the  unate  paradigm.  Typical  PLA's 
produce  shallow  recursion  trees  terminating  quickly  at  unate  leaves.  Espresso  II  uses 
this  benefit  to  trim  the  recursion  tree  and  balance  it  judiciously.  Besides,  using  matrix 
representation  for  representing  logic  function  and  minimization  process  in  Espresso  II 
requires  less  memory,  and  the  operations  to  be  executed  faster  than  other  forms  of  rep¬ 
resentation  of  function  [Ref.  5:  pp.  44-46]. 

'1  he  objectives  of  Espresso  II  are  to  minimize: 

•  The  number  of  the  product  terms  in  the  cover  (NPT) 

•  The  number  of  literals  (not  don't  care  )  in  the  input  parts  of  the  cover  (NLI) 

•  The  number  of  literals  in  the  output  part  (NLO) 

The  Espresso  II  minimization  procedure  defines  a  vector  objective  function 


<1>  =  (AT  7',  A  ’Ll,  MO) 


and  continues  to  iterate  through  into  main  minimization  loop  until  none  of  the  three 
components  of  O  have  been  reduced  since  the  last  past  through  the  loop  [Ref.  5:  pp. 
54-55] 
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Espresso  II  minimization  involves  seven  basic  routines  and  the  sequence  of  oper¬ 
ations  carried  out  by  Espresso  II  is  outlined  below: 

•  Complement :  Computes  the  complement  of  PLA's  and  the  don't  care  set. 

•  Expand  :  Expand  each  implicant  into  a  prime  and  remove  covered  implicant. 

•  Essential-Primes :  Extract  the  essential  primes  and  put  them  in  the  don't  care  set. 

•  Irredundant  Cover :  Find  a  minimal  (optionally  minimum)  irredundant  :over. 

•  Reduce  :  Reduce  each  implicant  to  a  minimum  essential  implicant. 

•  Iterate  :  Expand,  irredundant  cover,  reduce  until  no  improvement. 

•  Lastgrap  :  Try  reduce,  expand  and  irredundant  cover  one  last  time  using  a  different 
strategy.  If  succesful,  continue  the  iteration. 

•  Makesparse  :  Include  the  essential  primes  back  into  the  cover  and  make  the  PLA 
structure  as  sparse  as  possible. 

The  widely  description  and  explanation  of  over-all  Espresso  II  program  can  be  ob¬ 
tained  at  Ref  5, 
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APPENDIX  D.  AVERAGE  NUMBER  OF  PRODUCT  TERMS 


Table  3.  AVERAGE  NUMBER  OF  PRODUCT  TERMS 
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APPENDIX  E.  AVERAGE  COMPUTATION  TIME 


Table  4.  AVERAGE  COMPUTATION  TIME  (SEC) 


Number  of  Vari¬ 
ables 

Espresso  II 

MIH 

CISH 

1 

0.000 

0.000 

0.000 

2 

0.400 

0.890 

0.420 

3 

0.490 

5.760 

1.070 

4 

0.960 

6.010 

1.300 

5 

1.060 

6.480 

1.952 

6 

1.070 

6.860 

1.780 

7 

1.100 

7.976 

3.690 

S 

1.446 

11.410 

8.830 

9 

1.824 

24.704 

19.992 
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APPENDIX  F.  AVERAGE  MEMORY  USAGE 


Table  5.  AVERAGE  TOTAL  MEMORY  USAGE  (KBYTE) 


Number  of  Vari¬ 
ables 


Espresso  II 

MIH 

CISH 

3.600 

1.180 

2.160 

13.90 

7.660 

8.980 

64.160 

10.950. 

16.150 

77.000 

36.720 

45.786 

96.960 

57.980 

77.400 

103.200 

73.170 

100.22 

112.552 

115.090 

92.670 

127.080 

189.670 

137.060 

138.576 

319.696 

219.480 

APPENDIX  G.  PROGRAM  LISTING 


//include  <stdio.  h> 

//define  tnaxint  32767 

//define  maxvars  10  /*Max  //  of  variables  in  a  product  term.*; 

//define  maxcint  65535  /*2**( maxVars )-l*/ 

//define  maxcubes  1000  /*Max  //  of  cubes  to  allocate  per  level*/ 

//define  TRUE  1 

//define  FALSE  0 

//define  Function 

//define  max  5 

typedef  struct  cube  ] 

int  tt[ maxvars] ;  /*Bits  1  for  uncomplemented  variables.*/ 

int  ff[ maxvars];  /*Bits  1  for  complemented  variables.*/ 

int  ics;  /*  degree  of  clustering  of  minterm  */ 

int  oldies;  /*  degree  of  clustering  of  minterm  */ 

int  ic;  /*  minterm  coverage  of  each  implicant  */ 

int  selected;  /*  flag  of  selected  implicants  in  SOP  */ 

int  track;  /*  to  keep  track  of  implicants, cubes  */ 

jCube; 

int  numvar;  /*  number  of  the  actual  variables  */ 

int  m,n;  /*  counters  for  loops  */ 

int  aux.auxl;  /*  variables  for  temporary  impicants  */ 

int  naux.nauxl;  /*  variables  for  temporary  cubes  */ 

int  num_minterm;  /*  number  of  the  minterms  */ 

int  num_product;  /*  number  of  the  sum_of_product  terms  */ 

int  countimpl;  /*  counter  for  implicants  cover  cube  */ 

int  numcubes[ maxvars] ;  /*  flag  to  indicate  levels  and  cubes  */ 

int  covered[ maxvars] [ maxcubes] ;  /*  flag  to  indicate  covered  minterms  */ 

int  j,  k,kl,p,rc,r,rm,rl,rk,rz;  /*  Index  into  the  cubes .covered  array  */ 

int  found;  /*  flag  to  keep  track  of  covered  imp  */ 

int  isonum;  /*  counter  for  remainder  minterms  */ 

int  temp;  /*  dummy  variable  for  implicant  cube  */ 

int  small, high;  /*  variable  for  selection  sorting  */ 

int  num_ implicant;  /*  number  of  all  implicants  in  funct.  */ 

int  sub.subl;  /*  temporary  variables  */ 

int  dum.duml;  /*  temporary  variables  */ 

Cube  cubes] maxvars+1] [ maxcubes+l] ;  /*  Cube  representation  for  minterms  */ 

Cube  temporarycube;  /*  temporary  cube  for  manipulation  */ 

Cube  tempeube;  /*  temporary  cube  for  manipulation  */ 

Cube  implcubes] maxvars+1] [ maxcubes+1] ;  /*  cube  representation  for  impl.  */ 
Cube  tempi]  1]  [  max] ';  /*  temporary  implicant  */ 

Cube  nebor] 1] [ max] ;  /*  temporary  neighbor  cubes  */ 

FILE  *fpil,*fpol;  /*  pointers  for  the  files  */ 

Function  main  (arge,  argv) 
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int  argc; 
char  **argv; 

{ 

int  i,k; 
found  =  FALSE; 

for  (  m=0;  m  <=  maxvars  ;  m++  ) 

numcubesfm]  -  0;  /* 

for  (  ra=0;  m  <  numvar  ;  m++  ) 

for  (  i=0;  i  <  maxcubes  ;  i++  ) 

{ 


initialization  for  all  level 


covered! “Oil]  *  FALSE  ; 

cubes[ mj [ i] .  ics«0; 

cubesf m] t ij • oldics=0; 

cubes[  mj  [  ij .  ic=0; 

cubesf  mj  [  ij  •  selected  =0; 

cubes[  mj  [  ij .  track  =0; 

for  (  k=0;  k  <  numvar  ;  k++  ) 

cubes  [m]  [  i] .  tt[k]  *0; 

cubes!  m]  l  i]  •  f  f  [  k]  *»0; 


*/ 


/*  Read  the  minterm  from  input  file  */ 

readfile(argv[  1] ); 

for  (  m  »  0;  m  <  numvar;  m++  )  /*  For  all  level  except  the  last  */ 

for(j  *0;  j  <  numcubes[ m] ;  j44)  /*  For  all  cubes  at  this  level  */ 

for  (k  =j+l;  k  <  numcubes[m]  ;  k++)  /*  other  cubes  at  this  level*/ 

if (  rc=combinable(&cubes[ m] [ j]  ,&cubes[  m]  [  k]  ) ) 

covered! m] ! j]  =  TRUE;  /*  mark  the  cubes  as  covered  */ 

covered! m] [ k]  =  TRUE;  /*  mark  the  cubes  as  covered  */ 

/*  Combine  into  an  (m+l)-cube  store  in  tempcube  */ 

combine! &cubes!m]  [  j]  ,£tcubes[m]  [k]  ,&tempcube); 

found  =  FALSE  ;  /*  See  if  it  is  generated  before  */ 

for  (p=0;  p  <  numcubesfm+1] ;  p++) 

if  (  r  =  equalcubes(&cubes[  m+1]  [  p]  ,&tempcube)) 

found  =  TRUE; 

} 

) 

if( !  found) 

{ 
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/*  add  this  as  new  implicant  of  higher  level  */ 
for(kl=0;  kl  <numvar;  kl++) 

( 

cubes[  m+l]  [  numcubes[  m+1]  ] .  tt[  kl]  =  tempcube.  tt[  kl] ; 
cubes]  m+l]  [numcubes] m+l]  j .  ff[klj=  tempcube.  f f [ kl] ; 

} 

numcubes] m+l]  =numcubes[ m+l]  +  1; 


} 

} 

} 


/*  Find  all  possible  implicant  of  function  */ 
/* .  */ 

for  (m=0  ;m  <  numvar  ;  m++) 

{ 

for  (  j=  Oj  j  <  numcubes] m];  j++  ) 

{ 

if((m  ■*  0)  &&  ( !  (covered[  0]  [  j]  ))) 

{ 

num_impl leant ++; 
t  emp=num_imp 1 leant - 1 ; 

change(&cubes[m] [ j]  ,&implcubes[ 0] [ temp] ); 
if((m  1=  0)  &&  ( covered] m] [ j] )) 

l 

num_implicant++; 
temp=num_implicant - 1; 

change(&cubes[  m]  [  j] ,&implcubes[  0]  [  temp]  ); 

} 

if((m  !=0)  &&  (! covered] m] [ j] )) 

{ 

num_implicant++; 

temp=num_implicant-l; 

change ( Sicubes ]  m]  ]  j]  ,&implcubes[  0]  [  temp]  ); 

} 

1 


/*  find  the  Implicant  Cover  Size  (ICS)  */ 

/* . */ 


for(  tn=0;  m  <  num  minterm;  m++) 

{ 

for(  n=0;  n  <  num_implicant;'  n++) 

{ 

if(rm  =  imp_cov_size(&cubes[ 0] [ m] ,&implcubes[  0]  ] n]  ,&temporarycube)) 

1 

cubes]  Oj  [  m]  .  ics++;' 
cubes] 0] [ m] K oldics++; 
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} 


) 


/*  find  Minterm  Coverage  MC(I)  */ 

/* . */ 

for(m=0;  m  <  num_implicant  ;m++) 

[ 

implcubes] 0] [ m] . ic  =  coverage_size(&implcubes[ 0] [m] ); 


for(i=0;  i<  num_minterm; i++) 

i 

small=i; 

for( j=i+l; j<num  minterm; j++) 

{ 

if(cubes[ 0] [ j] .  ics<cubes[  0)  [ small] .  ics) { 
small^j; 

)} 

swap(&cubes[ 0]  [ small] ,&cubes[  0]  [ i]  ,&temporarycube); 

} 


/*  Application  of  CIS  to  find  minimal  sum  expression  */ 

/*  . */ 

fpol=fopen("cis.  o",'V'); 


REC:  for( i-0; i<num_m interm; i++) 

( 

if(cubes[ 0] [ 0] . ics  =  maxcint  &&  cubes[ 0] [ i] . ics  ! =  maxcint) 

{ 

swap(&cubes[ 0] [ i]  ,&cubes[  0]  [  0] ,&temporarycube); 
break; 


) 

for(i=0;  i  <  num  implicant;  i++) 

l 

high=i; 

for ( j=i+l;  j<num  implicant; j++) 

l 


if( implcubes[ 0] [ j] . ic  >  implcubes] 0] [ high] .  ic) { 

swap(&implcubes[ 0] ( j] ,&implcubes[ 0]  [high] ,&temporarycube); 

n 

} 


/*  Extended  Search  Technique  */ 

/* . */ 


if((cubes{ 0]  [  0] . ics  =  2)) 

[ 

aux=0;' 

auxl=0; 
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sub=0; 

subl=0; 

naux=0; 

nauxl=0; 

for(m=0; m<num_impl leant; m++) 

[ 

i£(rl=cover(&cubes[  0]  [0]  ,&implcubes[ 0]  [in]  ,&tempcube)) 
sub++; 

dunFs  ub  -  1; 

swap(&tempi[  0]  [  dum]  ,&implcubes[ 0] [m] ,&tempcube); 
tempi[ 0] [ dum] .  track  =  m; 


/* 


if(teropi[  0)  [  0] .  ic  >  tempi[  0]  [  1] .  ic) 

{ 


tempi[  0]  [  0] .  selected  *  1; 
aux*tempi[  0]  [  0] .  track; 


auxl*tempi[ 0] 
swap(&tempi[ 0 
swap(&tempi[  0 


1] .  track; 

[ 0] ,Sdmplcubes[ 0] [ aux] ,&temporarycube); 

[ 1] ,&implcubes[ 0] [ auxl] ,&temporarycube); 


*/ 


for( j=0; j<num  minterm; j++) 

[ 

if((rz«cover(&cubes[ 0] [ j] ,&implcubes[ 0] [ aux] ,&temporarycube))  && 
(cubes[  0] [ j] . ics  !=  maxcint  )) 

cubes[ 0] [ j] . ics  »  maxcint; 
isonum++; 

for(k*0;  k<num_implicant; k++) 

{ 

if ( ( r=cover(&cubes[  0]  [  j]  ,&implcubes[ 0] [ k] ,&temporarycube) )  && 
( implcubesf 0] [ kj .  selected  !=  1)  &&  ( implcubes[ 0] [ k] . ic>0)) 

--implcubes[ 0] [ k] . ic; 

} 

} 

) 

) 

if(isonum  ==  num  minterm) 

{ 

goto  EXIT; 

1 

goto  REC; 


/* 


*/ 

i 


if(tempi[  0]  [  1]  .  ic  >  tempi[  0]  [  0]  .  ic) 

[ 

tempi [  0]  [ 1] .  selected  =  1; 
aux=tcmpi[  0]  [  1]  .  track; 
auxl=tempi[ 0]  (  0] , track; 

swap(&tetr.pi(  0]  [  1]  ,&implcubes[  Oj  [  aux]  ,&tempcube); 
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swap(&tempi[ 0] [ 0]  ,&implcubes[ 0]  [ auxl] ,&tempcube); 
for(  j=0;  j<num_minterm;  j++) 

if((rz=cover(&cubes[  0]  {  j]  ,& imp 1 cubes [ 0] [ aux] ,&tempcube))  && 
(cubes] 0]  [  j] .  ics  !=  maxcint  )) 

{ 

cubes]  0}  [  ,i) .  ics  =  maxcint; 
isonum++; 

for(k=0;  k<num_implicant;  k++) 

if((r=cover(&cubes]  0]  [  j] ,&implcubes[ 0] [k] ,&tempcube))  && 
(implcubes]  0]  ]  k] .  selected  !=  1)  &&  (implcubes]  0]  ]k] .  ic  >0)) 

{ 

--implcubes]  0]  [  k] .  ic; 

} 

} 

) 

i 

if(isonum  =  num_rainterm) 

{ 

goto  EXIT; 

} 

goto  REC; 

} 


/* 


*/ 


if ( (tempi[  0]  [  0] .  ic  «  1)  &&  (tempi!  0]  [  1] .  ic  «*  1)) 

! 

for(m=0; m<num_minterm;  m++) 

if(r  l=combinable(&cubes[  0]  [  ni]  ,&cubes] 0] [ 0]  ,&tempcube)) 

l 

subl++; 

duml=subl  -  1; 

swap(&nebor] 0] ( dural]  ,&cubes[ 0]  [  m]  ,&tempcube); 
nebor{  0]  [  dural] .  track  =  m; 

1 

1  /* . */ 


if((nebor[ 0] [ 0] . oldies  =  2)&&(nebor] 0] [ 1] . oldies  !=  2)) 

! 

if(rl=cover(&nebor[  0]  f 0]  ,&tempi] 0] ] 0] ,&tempcube)) 

! 

tempi] 0] [ 0] . selected  =  1; 
aux=terapi( 0] { 0] . track; 
auxl=tempi[  0]  (  1] .  track; 

swap(&tempi[ 0] [ 0] j&implcubes] 0] ] aux] ,&temporarycube); 
swap(&tempi[ 0]  |  1]  .Simplcubss] 0] [ auxl]  ,&temporarycube); 
naux=nebor( 0] ( 0] . track; 
nauxl=nebor( 0] ] 1] .  track; 

swep(6nebor[ 0] ( 0] ,&cubes[ 0] { naux] ,&temporarycube); 
swap(&nebor( 0]  ]  1]  ,&.cubes[0]  [nauxl]  ,&temporarycube); 
for(  j=0;  j<nuir._m  in  terra;  j++) 
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if((rz=cover(&cubes[ 0]  [ j] ,&implcubes[ 0] [ aux]  ,&temporarycube))  && 
(cubes] 0] [ j] . ics  !=  maxcint  )) 

{ 

cubes[ 0] [ j] . ics  =  maxcint; 
isonum-H-; 

for(k=0;  k<num_impl leant; k++) 

{ 

if((r=cover(&cubes[ 0] [ j] ,&implcubes[ 0] [k] ,&temporarycube))  && 
(implcubes] 0]  [k] .  selected  !=  1)  &&  ( imp 1 cubes [ 0]  [  k] .  ic  >  0  )) 

{ 

--implcubes]  0]  [k] .  ic; 

) 

} 

} 

} 

if(isonun)  —  num_minterm) 

{ 

goto  EXIT; 

} 

goto  REC; 

} 

else 

l 

tempi[  0]  [  1] .  selected  *  1; 
aux=tempi[  0]  [  1] .  track; 
auxl=tempi[ 0] [ 0] .  track; 

swap(&tempi[  0]  [  1]  ,&implcubes[  0]  [  aux]  ,6ftempcube); 
swap(&tempi[  0]  [  0] ,&implcubes[ 0] [ auxl] ,&terapcube); 

nau.\=nebor[  0]  [  0] .  track; 
nauxl=nebor[ 0]  [  1] .  track; 

swap(&nebor[  0]  [  0] ,&cubes[ 0J [ naux] ,&temporarycube); 
swap(&nebor[  0]  (  1]  ,6(Cubes[  0]  [  nauxl]  ,&temporarycube); 
for( j=0; j<num  minterm; j++) 

{ 

if((rz=cover(&cubes[ 0] [ j] ,&implcubes[ 0] [ aux] ,&tempcube))  && 
(cubes[ 0] [ j] , ics  !=  maxcint  )) 

{ 

cubes[ 0] [ j] . ics  =  maxcint; 

' sonum++; 

>or(k=0; k<num  implicant; k++) 

( 

if(( r=cover(&cubes[ 0] [ j] ,&implcubes[ 0] [ k] ,&tempcube))  && 

( implcubes[  0]  l k] .  selected  !=  1)  &&( imp 1 cubes [ 0] [ k] . ic>0)) 

{ 

--implcubes[ 0] [ k] . ic; 

! 

} 

) 

i 

if(isonum  ==  num  minterm) 
goto  EXIT; 

] 

goto  REC; 
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if((nebor[ 0}  [ 1] . oldies  —  2)  &&  (nebor[ 0]  [  0] .  oldies  !=  2)) 

{ 

if ( r l=cover(&nebor[ 0] [ 1] ,&tempi[ 0] [ 1] ,&tempcube) ) 

[ 

tempif  0]  [  1] .  selected  =  1; 
aux=tempi[ 0]  [ 1] .  track; 
auxl=tempi[  0]  [  0] .  track; 

swap(&tempi[ 0] [ 1] ,&implcubes[ 0] [ aux] ,&tempcube); 
swap(&tempi[ 0]  [ 0] ,&implcubes[ 0] [ auxl] ,&tempcube); 

naux=nebor[ 0]  [  0] . track; 
nauxl=nebor[  0J [ 1] . track; 

swap(&nebor| 0]  [ 0] ,&cubes[ 0] [ naux] ,&temporarycube); 
swap(&nebor[  0]  [  1]  ,6icubes[  0]  [  nauxl]  ,&temporarycube); 
for( j=0; j<num_m interm; j++) 

{ 

ifUrz^coverC&cubesf  0]  [  j] ,& imp 1 cubes [ 0]  [ aux] ,&tempcube))  && 

( cubes [ 0] [ j] . ics  !  =  maxcint  )) 

[ 

cubes[ 0] [ j] . ics  =  maxcint; 
isonum++; 

for(k=0; k<num_implicant; k++) 

if((r=cover(&cubes[ 0] [ j] ,&implcubes[ 0]  [  k] ,&tempcube))  && 

( implcubes[ 0] [ k] . selected  Is  1)  &&  (implcubes[ 0] [k] . ic>0)) 

{ 

--implcubes[ 0] [ k] . ic; 


} 

} 

if(isonum  =  num  minterm) 

{ 

goto  EXIT; 

) 

goto  REG; 

) 

else 

{ 

tempi(  0]  [  0J .  selected  =  1; 
aux=terapi(  0]  [  0] .  track; 
auxl=tempi[ 0]  (  1]  .  track; 

swap(Sctempi[  0]  [  0J  ,&implcubes[  0]  [  aux]  ,&temporarycube); 
swap(&tempi[ 0]  [ 1]  ,&implcubes[ 0] [ auxl]  ,&temporarycubej; 

naux=nebor[  0]  [  0] .  track; 
nauxl=nebor[  0]  [  1]  .  track; 

swap(&nebor| 0]  [0]  ,&cubes[0] [naux] ,&temporarycube); 
swap(&nebor[ 0] [ 1] ,&cubes[ 0] [ nauxl] ,&temporarycube); 
for(  j=0;  j<num_minterin;  j++) 


if((rz=cover(&cubes[ 0] [ j] ,&implcubes[ 0] [ auxj  ,&temporarycube))  && 
(cubes] 0]  [  j] .  ics  J =  maxcint  )) 

{ 

cubes[0] [ j] . ics  =  maxcint; 
isonum++; 

for(k=0; k<num_implicant;  k++) 

{ 

if((r=cover(&cubes[ 0] [ j]  ,&implcubes[ 0]  [  k]  ,&temporarycube))  && 
(implcubes] 0] [k] .  selected  !=  1)  &&  ( implcubes] 0]  [  k] .  ic>0)) 

[ 

--implcubes]  0]  [k] .  ic; 

} 

} 

) 

} 

if(isonum  “  num  minterm) 

{ 

goto  EXIT; 

J 

goto  REC; 

} 

} 


/* . */ 

if((nebor[ 0] [ 0] . oldies  ! =2  &&  nebor[ 0] [ 1] . oldies  ! =2) 

(nebor[  0]  l  0] .  oldies  =2  &&  nebor]  0]  [  1] .  oldies  =2)) 

{ 

tempi[  0]  [  0] .  selected  =  1; 
aux=tempi]  0]  [  0] . track; 
auxl^tempi]  0]  [  1] .  track; 

swap(&tempi[  0]  ]  0]  ,&implcubes[  0)  [  aux]  ,&temporarycube); 
swap(&tempi]  0]  [  1]  ,&implcubes]  0]  ]  auxl]  ,&tempor  ary  cube); 
naux=nebor[ 0] ( 0] . track; 
nauxl=nebor[ 0] [ 1] . track; 

swap(&nebor[  0]  [ 0]  ,&cubes( 0] [ naux] ,&temporarycube); 
swap(&nebor[ 0] [ 1] ,&cubes[ 0] [ nauxl] ,&temporarycube); 
for( j”0;  j<num_minterm; j++) 

{ 

if( ( rz=cover(&cubes] 0] [ j] mpl cubes [ 0] [ aux] ,&temporarycube))  && 
(cubes]  0]  ]  j] ics  !=  maxcint;  )) 

{ 

cubes] 0] [ j] . ics  =  maxcint; 
isonu:n++; 

for(k=0; k<num  implicant; k++) 

{ 

if((r=cover(&cubes] 0] { j] ,& imp 1 cubes] 0J [ k] ,&temporarycube))  && 
( implcubes] 0] ] k] , selected  !=  1)  &&  ( implcubes] 0]  ( k] .  ic>0)) 

{ 

--implcubes] 0  ( k] . ic; 

J 

I 

) 

} 

if(isonum  ==  num  minterm) 

{ 
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goto  EXIT; 

} 

goto  REC; 

} 

} 


/* 


*/ 


if(( tempi]  0)  [  0] .  ic  —  2)  &&  (tempi[  0]  [  1] .  ic  =  2)) 

{ 

for(m=0;  m<num_minterm;  m++) 

{ 

if ( rl=combinable(&cubes[ 0] ] 0] ,&cubes[ 0] [ m] ,&terapcube) ) 

l 

subl++; 

duml=subl  -  1; 

swap(&nebor[  0]  [duml]  ,&cubes[0]  [m]  ,&tempcube); 
nebor[  0]  [  duml] .  track  =  m; 

} 

1 

if((nebor[ 0]  [  0] .  oldies  —  2)  &&  (nebor[ 0] [ 1] . oldies  !=2)) 

{ 

if(rl=cover(&nebor[ 0] [ 0]  ,&tempi[  0]  [  0]  ,&tempcube)) 

{ 

tempi[  0]  [  0] .  selected  =  1; 
aux=tempi[ 0] J  0] . track; 
auxl=tempi[  0] [ 1] .  track; 

swap(&tempi[ 0] [ 0] ,&implcubes[ 0] [ aux] ,&temporarycube); 
swap(&tempi[  0] [ 1] ,&implcubes[ 0] [ auxl] ,&temporarycube); 

naux=nebor[ 0] [ 0] . track; 
nauxl=nebor[ 0] [ 1] . track; 

swap(&nebor[ 0] [ 0] ,&cubes[ 0] [ naux] ,&tempor ary cube); 
swap(&nebor[ 0] [ 1] ,&cubes[oj [nauxl] ,&temporarycube); 

for( j=0; j<num_minterm;  j++) 

{ 

if ( (rz=cover(&cubes[ 0] [ j] ,& imp 1 cubes [  0]  [  aux] ,&temporarycube))  && 

( cubes] 0] [ j] . ics  !=  maxcint  )) 

{ 

cubes] 0] [ j] . ics  =  maxcint; 
isonutn-H-; 

for(k=0; k<num_implicant;  k++) 

i 

if ( ( r=cover(&cubes[ 0] ] j] ,&implcubes[ 0]  [  k] ,&temporarycube))  && 
( implcubes] 0] [ k] . selected  !=  1)  &&  ( implcubes] 0]  [ k] .  ic>0) ) 

{ 

-•implcubes] 0] ] k] . ic; 

} 


1 

! 

if(isonum  —  nurn  minterm) 

{ 

goto  EXIT; 
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} 

goto  REC; 

} 

else 

{ 

tempi[  0]  [  1] .  selected  =  1; 
aux=tempi[  0]  [  1] .  track; 
auxl=tempi[ 0] [ 0] .  track; 

swap(&tempi[  0]  [  1]  ,&implcubes[  0]  [  aux]  ,&tempcube); 
swap(&tempi[  0]  [  0]  ,&implcubes[  0]  [ auxl]  ,&tempcube); 

naux=nebor[ 0]  (  0] .  track; 
nauxl=nebor[ 0]  [ 1] .  track; 

swap(&nebor[ 0] ( 0] ,&cubes[ 0] [ naux]  ,&temporarycube); 
swap(&nebor{ 0] [ 1]  ,&cubes[  0]  [nauxl]  ,&temporarycube); 
for( j=0; j<num_rainterm; j++) 

{ 

if((rz=cover(&cubes[  0]  {  j]  ,&implcubes[  0]  [  aux]  ,&tempcube))  && 
(cubes] 0] [ j] . ics  !  =  maxcint  )) 

[ 

cubes! 0]  [ j] .  ics  -  maxcint; 
isonum++; 

for(k=0;  k<num  implicant;  k++) 

{ 

if((r-cover(&cubes[ 0] [ j]  ,&implcubesf 0]  [  k]  ,&tempcube))  && 

( implcubesj 0]  [  k] .  selected  !=  1)  &&  (implcubes[  0]  [k] .  ic>0)) 

! 

--iraplcubes] 0] [ k] . ic; 

} 

} 

} 

) 

if(isonum  =  num_m interm) 

! 

goto  EXIT; 

} 

goto  REC; 

) 

} 


/* . */ 

if((nebor[ 0] [ 1] . oldies  =  2)  &&  (nebor! 0]  ( 0] . oldies  !=2)) 

{ 

if ( rl=cover(&nebor[ 0] [ 1] ,&tempi[ 0] ! 0] ,&tempcube) ) 

{ 

tempi[  0]  [  0]  ,  selected  =  1; 
aux=tempi[ 0]  [  0]  .  track; 
auxl=tempi! 0]  [  1]  .  track; 

swap(&tempi[ 0] [ 0] ,&implcubes[ 0] [ aux]  ,&temporarycube); 
swap(&tempi{  0]  [  1]  ,&implcubesj  0]  [  auxl]  ,&temporarycube); 

naux=nebor[ 0]  {  0]  .  track; 
nauxl=nebor!  0]  [  1] ...  track; 

swap(&nebor[  0]  [  0]  ,5tcubes[  0]  [  naux]  ,&tempofarycube); 
swap(&nebor[  0]  [  1]  ,&cubes[  0]  ( nauxl]  ,&temporarycube); 
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for( j=0; j<num_minterm; j++) 

l 

if((rz=cover(&cubes[ 0]  [  j]  ,&implcubes[ 0] [ aux] ,&temporarycube))  && 
(cubes] 0]  [  j] .  ics  !=  maxcint  )) 

{ 

cubes[0] [ j] .  ics  =  maxcint; 
isonum++; 

for(k=0;  k<num  imp li cant ; k++) 

{ 

if((r=cover(&cubes] 0] [ j] ,&implcubes]  0]  [k] ,&temporarycube))  && 
(implcubes] 0]  [ k] .  selected  !=  1)  &&  ( implcubes] 0]  [k] .  ic>0)) 

{ 

--implcubes] 0] [ k] . ic; 

} 

} 

} 

} 

if(isonum  —  num  minterm) 

{ 

goto  EXIT; 

} 

goto  REC; 

} 

else 

{ 

tempi[  0]  [  1] .  selected  =  1; 
aux=tempi[  0]  [  1] .  track; 
auxl=tempi[ 0] [ 0] . track; 

swap(&tetnpi[  0]  [  1]  ,&implcubes[  0]  [  aux]  ,&tempcube); 
swap(&tempi[ 0] [ 0] ,&implcubes[ 0] [ auxl] ,&tempcube); 

naux=nebor[ 0] [ 0] . track; 
nauxl=nebor[  0J  [  1] .  track; 

swap(&nebor[  0]  [0]  ,£ccubes[0]  [naux]  ,&temporarycube); 
swap(£tnebor[  0]  [  1]  ,&cubes[oj  [nauxl]  ,&temporarycube); 
for( j=0; j<num_ra interm;  j++) 

{ 

if ( ( rz=cover(&cubes[ 0] [ j] ,&implcubes[ 0] [ aux] ,&tempcube) )  && 
(cubes] 0] [ j] . ics  !=  maxcint  )) 

{ 

cubes] 0] [ j] . ics  =  maxcint; 
isonum++; 

for(k=0; k<num_implicant;  k++) 

{ 

if ( ( r=cover(&cubes[ 0] [ j] ,&implcubes] 0] [ k] ,&tempcube))  && 

( implcubes] 0] ] k] . selected  !=  1)  &&  ( implcubes] 0] [ k] .  ic>0)) 

{ 

--implcubes] 0] ( k] . ic; 

1 

} 


] 

if(isonum  —  num_minterm) 

{ 


goto  EXIT; 

} 
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goto  REC; 


} 

} 

/* . */ 

if((nebor[  0]  [  0] .  oldies  —  2  &&  nebor[  0]  [  1) .  oldies  =2) 

(nebor[  0]  [  0] .  oldies  !=  2  &&  neborf  0]  [  1] .  oldies  !=2)) 

{ 

tempif  0]  [  0] .  selected  =  1; 
aux=tempi[ 0]  [  0] .  track; 
auxl=tempi[  0]  [  1] .  track; 

swap(&tempi[  0]  [  0]  ,&itnplcubes[  0]  [  aux]  ,&temporarycube); 
swap(&tempi[ 0] [ 1] ,&implcubes[ 0]  [ auxl] ,&temporarycube); 

naux=nebor[  0]  [ 0] .  track; 
nauxl=nebor[  0]  [  1] .  track; 

swap(&nebor[ 0]  [  0]  ,&cubes[  0] { naux] ,&temporarycube); 
swap(&nebor[ 0] [ ]] ,&cubes[0] [nauxl] ,&temporarycube); 

£or(  j=0;  j<nutn  niinterm;  j++) 

{ 

if((rz=cover(&cubes[ 0J [ j] ,& imp 1 cubes [ 0] [ aux] ,&temporarycube))  && 
(cubes[ 0] [ j] .  ics  !=  maxcint  )) 

[ 

cubes[ 0]  [ j] .  ics  «  maxcint; 
isonum++; 

£or(k=0; k<num_implicant; k++) 

[ 

if((r=cover(5tcubes{  0]  [  j]  ,&implcubes[  0]  [  k]  ,&temporarycube))  && 
( implcubesf 0] [ k] . selected  !=  1)  &&  ( implcubes[ 0]  [ k] .  ic>0)) 

{ 

-”implcubes[  0]  [  k]  ■<  ic; 

} 

} 

} 

} 

if(isonum  ==  nutn  minterm) 

[ 

goto  EXIT; 

1 

goto  REC; 

} 

) 

) 


/* 


*/ 


if(cubes[ 0] [ 0] . ics  !  =2) 

I 

for(m=0;  m<num_implicant; m++) 

{ 

if( rc=cover(&cubes[  0]  [  0]  ,5iiinplcubes[  0]  [  m]  ,&temporarycube)  && 

(  implcubesf  0]  [  m] ...  selected  !=  1)  &&  (  implcubes[  0]  [  m] ic  >0)) 
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[ 

implcubes[ 0] [m] .  selected  =  1; 
implcubesj 0]  [in] .  ic  =  0; 
for( j=0; j<nuro_minterm;  j++) 

{ 

if ( rz=cover(&cubes[ 0] [ j] ,&implcubes[ 0] [ m] ,&temporarycube) 

&&  (cubes[ 0]  [  j] .  ics  !=  maxcint)) 

{ 

cubes[ 0]  [  j] .  ics  =  maxcint; 
isonum++; 

for(k=0;  k<num_itnplicant;  k++) 

{ 

if(r=cover(&cubes[ 0] [ j]  ,&implcubes[  0]  [  k]  ,&temporarycube) 
&&  (implcubes[  0] [k] .  selected  !=  1)&& 

( implcubes[ 0]  [ kj . ic  >  0  )) 
f 

--implcubes[ 0] [ k] . ic; 

} 

) 

} 

} 

if(isonum  =  num  minterm) 

{ 

goto  EXIT; 

} 

goto  REC; 

) 

} 


/*  - . . . . . .  */ 

EXIT:  for(m=0; m<num_implicant; m++) { 

if C implcubes[  0]  [  m]  .  selected  =  1) 

{ 

printcube(&implcubes[ 0) [ m] ); 
num  product++; 

} 

} 

fclose( fpol); 

printf("%d  n",num  product); 

) 


/* . */ 


Function  imp_cov_size(al,a2,a3) 
Cube  *al,*a2J*a3; 

[ 

int  i; 
int  v_ics; 
int  check; 
check  =0; 
v_ics  =0; 
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/*  check  tnat  this  minterm  can  be  covered  by  this  implicant  or  not  */ 
for(i  =0;  i  <  numvar  ;  i++)  [ 
a3  ->tt[  i]  =  AND(al->tt[  i]  ,  a2->tt[i]); 
a3  ->ff[i]  =  AND(al->ff[i]  ,  a2->ff[  ij ); 
if((a3->tt[i]  =  a2->tt{  i] )  &&  (a3->ff[i]  =  a2->ff[i])) 

check++; 

}} 

/*  if  all  bits  are  same  as  implicant  than  retuns  as  v_ics  =1  */ 

/*  Otherwise  it  remains  as  0.  */ 

if C  check  =  numvar) 

{ 

v_ics=l; 

} 

return(v_ics); 


/* 


*/ 


Function  int  coverage_size(al) 

Cube  *al; 

[ 

int  i;' 

int  parameter; 

int  messanger; 

parameter  =  0; 

for(i=0;  i  <  numvar;  !+■+){ 

i£((al->tt[i]  ==  0)  &&  (  al->ff[ i]  =  0)){ 

parameter++;  J  ] 

if(parameter  =  0){ 
messanger  “1; j 

if(parameter  ==  1){ 
messanger  =  2; j 

if(parameter  —  2){ 
messanger  =4; } 

if (parameter  ==  3){ 
messanger  =8;  j 

if(parameter  ==  4){ 
messanger  =  16;  j 

if(parameter  =  5){ 
messanger  =  32;  j 

if(parameter  — 6){ 
messanger  -  64; } 


69 


if(parameter  =  7){ 
messenger  =  128; } 

if(parameter  *=»  8){ 
messenger  =  256; } 

if(parameter  ==  9){ 
messenger  =  512; } 

return( messanger ) ; 


/* . 

Function  int  swap(al,a2,a3) 

Cube  *al,'va2,*a3; 

{ 

int  i; 

for(  is0;  i  <  numvar  ;  i++) 

[ 

a3->tt[i]  *  al->tt[  i] ; 
al->tt[ ij  =  a2->tt[  i] ; 
a2->tt|ij  «=a3->tt[ij; 

a3->ff[  i]  =*  al->ff(i); 
al->ff[i]  ■  a2->f £[  i] ; 
a2->ff[i]  -  a3->ff[i]; 

} 

a3->ics  =  al->ics; 
al->ics  =  a2->ics; 
a2->ics  =  a3->ics; 

a3->oldics  =  al->oldics; 
al->oldics  =  a2->oldics; 
a2*>oldics  =  a3->oldics; 

a3->ics  =  al->ics; 
al->ics  =  a2->ics; 
a2~>ics  =  a3->ics; 

a3->selected  =  al->selected; 
al->selected  =  a2->selected; 
a2->selected  =  a3->selected; 

a3->track  =  al->track; 
al->track  =  a2->track; 
a2->track  =  a3->track; 


} 


/* 


’0 


*/ 


Function  int  cover(al,a2,a3) 

Cube  *al,*a2,*a3; 

{ 

int  i; 

int  v_cover; 
int  check; 
check=0; 
v_cover=0; 

£or(i=0;  i  <  numvar;  i++){ 

a3->tt[ i]  =  AND(al->tt[ i] ,a2->tt[ i] ); 
a3->ff[  i]  =  AND( al->ff[  i]  ,a2->ff[  i] ); 
if ( ( a3->tt[  i]  =  a2->tt[  i] )  &&  (a3->ff[i]  —  a2->ff[i])) 
{ 

check++; 

1 

1 

if (  check  **  numvar) 

{ 

v_cover=l; 

} 

return(v  cover); 

} 


/* 


*/ 


Function  int  equalcubes(al,a2) 

Cube  *al,*a2; 

{ 

/*  if  EQ  return  1  else  return  0  */ 
int  i; 

int  v_equal; 
v_equal  =1; 

for  (i  =  0;  i  <numvar  ;  i++)  { 

if  ( ( a  1  ->tt[  i]  !  =  a2->tt[  ij  al->ff[i]  !  =  a2->ff[  i] )) 

{ 

v_equal  =  0; 

} 

! 

return(v  equal); 

) 


/* . */ 


Function  int  combinable(al,a2) 
Cube  *al,"a2; 


{ 

int  i; 

int  v_combinable; 
int  cone; 

cone=0;  /*  1:  combinable  else  not  */ 

v_combinable  =0; 

for  (i  =  0;  i  <numvar  ;  i++)  { 

if  C((al->tt[i]  !  =  a2->tt[  i] ) )  ((al->ff[i]  !=  a2->ff[  i]  ))) 

{ 

cone++; 

} 

} 

if  (cone  =1)  v_combinable=l; 
r eturn( v_combinab 1 e ) ; 

} 


/* 


*/ 


Function  reodfile( filename) 
char  ’’'filename; 

{ 

int  i, j; 

Cube  *cl;  /*  cl  point  to  c[0][0]  */ 

if((fpil=fopen(filename,"r"))  —  NULL) 

printf("  nERROR  -  Cannot  open  designated  read  file  n"); 
return; 

} 

while(! feof(fpil)){ 


fscanf( fpil ,"%d  %d  n",  &num_minterm,&numvar); 
numcubes[0]  =  num_minterm; 
for(i=0;i  <  num_minterm;  i++) j 
for(j=0; j  <  numvar;  j++) 

fscanf(fpil ,"%ld" ,&cubes[ 0]  [  i] .  tt[ j] ); complement(&cubes[  0]  [  i]  ); 

} 

) 

) 


/*  . */ 


Function  int  complement(al) 
Cube  *al; 

{ 

int  j; 

for(j=0;  j  <  numvar;  j++) 

{ 

if  (al  ->  tt[j]  ==  1) 

{ 
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al  ->  ff[j]  =  0; 

1 

else 

al  ->  ff[ jj  =1  ; 


/* . */ 


Function  cotnbine(al,a2,a3) 

Cube  ,val,*a2,*a3  ; 

{ 

int  i; 

for  (i  *  0;  i  <maxvars  ;  i++)  [ 

a3«>tt[i)  ■  AND(al->tt[i]  ,  «2->tt[i])  ; 

a3->ff[  i]  »  AND(al->ff[i]  ,  a2->ff[i]j  ; 

I 


/* 


*/ 


Function  int  AND(il,i2) 
int  il , 12; 

I 

if(il=*12)  return(il); 
else  return(O); 

) 


*/ 


Function  Ant  printcube( al) 

Cube  *al  ; 

1 

int  j; 

for( j=0; j  <  numvar;  j++) { fprintf( fpol ,"%ld",(al*>tt[ j) )); } 

fprintf( ipol n'); 

for( j=0; j<numvar; j++) { fprintf( fpol , al*>f f[  j]  );  j 

fprintff  fpol n"); 
fprintf( fpol n"); 

} 


/* 


*/ 
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Function  int  change(al,a2) 
Cube  *al,*a2; 

{ 

int  i; 

for(  i=0;  i<  numvar;  i++) 

{ 

a2->tt[  i]  =al->tt[i]; 
a2->ff[i]  =  al->ff[  i] ; 


74 


V 


LIST  OF  REFERENCES 


1.  P.  Tirumalai  and  J.  T.  Butler,  Minimization  Algorithms  for  Multiple-Valued  Pro¬ 
grammable  Logic  Arrays,  Forthcoming  in  IEEE  Transactions  on  Computers. 

2.  J.  K.  Breeding,  Digital  Design  Fundamentals,  Prentice  Hall,  Englewood  ClifFs,  New 
Jersey,  1989. 

3.  J.  F.  Wakcrly,  Digital  Design  Principles  and  Practices Prentice  Hall,  Englewood 
Clift's  ew  Jersey,  1990. 

4.  Y.  Igaraski,  An  Improved  Lower  Bound  on  the  Maximum  Number  of  Prime 
lmplicants ,  The  Transactions  on  the  IECE  of  Japan,  Vol  E62,  June  1979, 
pp. 389-394. 

5.  R.  K.  Bravton,  Gary  D.  Hacht,  Logic  Minimization  Algoritms  for  VLSI  Synthesis, 
Kluwer  Academic  Publishers,  19S4. 

6.  S,  H.  Unger,  The  Essence  of  Logic  Circuits,  Prentice  Hall,  Englewood  Clifts,  New 
Jersey,  1989. 

7.  U.  Manber,  Introduction  to  Algorithms,  Addison-Wesley,  Menlo  Park,  1989. 


73 


8.  P.  Tirumalai  and  J.  T.  Butler,  Analysis  of  Minimization  Algorithms  for  Multiple- 
Valued  Programmable  Logic  Arrays,  Proceeding  of  18th  International  Symposium 
on  Multiple- Valued  Logic,  May  1988  pp.  226-236. 

9.  Chyan  Yang  and  Y.  M.  Wang,  A  Neighborhood  Decoupling  Algorihm  for  Truncated 
Sum  Minimization,  Proceeding  of  the  20th  International  Symposium  on  Multiple- 
Valued  Logic,  May  1990, 

10.  E.  J.  McCluskcy,  Logic  Design  Principles  with  Emphasis  on  Testable  Semicustom 
Circuits,  Prentice  Hall,  Englewood  Cliffs,  New  Jersey,  1986. 

11.  H.  Taub,  Digital  Circuits  and  Processors,  McGraw-Hill,  NY,  1982, 


76 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5002 

3.  Deniz  Kuvvetleri  K.ligi  1 

Personel  Sube  Bsk.ligi 

Bakanliklar,  Ankara  /  TURKEY 

4.  Kara  Harp  Okulu  K.ligi  1 

Kutuphanesi 

Bakanliklar,  Ankara  /  TURKEY 

5.  Deniz  Harp  Okulu  K.ligi  1 

Kutuphanesi 

Tuzla,  Istanbul  /  TURKEY 

6.  Hava  Harp  Okuiu  K.ligi  1 

Kutuphanesi 

Ycsilyurt,  Istanbul  f  TURKEY 

7.  Orta  Dogu  Teknik  Universitesi  1 

Okul  Kutuphanesi 

Balgat,  Ankara  /  TURKEY 

8.  Bogazici  Universitesi  1 

Okul  Kutuphanesi 

Bebek,  Istanbul  l;  'I  URKLY 

9.  Ege  Dokuz  Evlul  Universitesi  1 

Okul  Kutuphanesi 

Cumhuriyet  .Vleydani,  Izmir TURKEY 

10.  Chairman,  Code  EC  1 

Department  of  Electrical  and  Computer  Engineering 

Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

1 1 .  Professor  Ch\  an  Yang  2 


Department  of  Electrical  and  Computer  Engineering,  Code  EC  Ya 
Nasal  Postgraduate  School 
Monterey,  CA  93943-5000 


Professor  Jon  T.  Butler 

Department  of  Electrical  and  Computer  Engineering,  Code  EC/Bu 
Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

Ugur  Ozkan 

Had  Zihni  Efendi  Sok.  No.  8/3 
Capa,  Istanbul  /  TURKEY 

Kadri  Hekimoglu 
DeGol  Caddesi  No.  1/5 
Tandogan,  Ankara  /  TURKEY 


